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Abstract. This final article about the CHOOZ experiment presents a complete description of the V c source 
and detector, the calibration methods and stability checks, the event reconstruction procedures and the 
Monte Carlo simulation. The data analysis, systematic effects and the methods used to reach our conclu- 
sions are fully discussed. Some new remarks are presented on the deduction of the confidence limits and 
on the correct treatment of systematic errors. 



1 Introduction 
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as follows: 
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Neutrino oscillation experiments are sensitive probes of 
the possible existence of a finite neutrino mass and pro- 
vide a way to study physics beyond the Standard Model 
of electroweak interactions 1:. In fact, lepton flavour vi- 
olation and the existence of nonzero neutrino masses can 
give rise to neutrino oscillations, as first pointed out by 
Pontecorvo and Maki et al. A_ . Several experiments, 
studying solar |5]|5][7||H] or atmospheric neutrinos |91ll0lllll 
I12U13II14) . have measured fluxes consistently lower than 
expectations. This can be interpreted as due to various 
forms of neutrino oscillations. In particular the so-called 
"atmospheric neutrino anomaly" is the observation of a 
Vfj,/v e ratio which is roughly one half of what expected 
and its possible explanation might be due to either os- 
cillation of «-> v T or to v ^ <-> v c . In a model with two 
neutrino eigenstates of mass m\ and 7712 which mix to form 
two flavour states, a pure beam of electron-flavoured neu- 
trinos has a survival probability which oscillates due to 
the mi — TO2 mass difference. For a single neutrino energy 
E v (MeV) and a distance from the source L (meters), the 
survival probability can be written in terms of the mixing 
parameter sin 2 28 and the difference of the squared masses 



P(F e -> V e ) = 1 - sin 2 26» sin 2 



1.27<STO 2 (eV 2 )L(m) 
E„{MeV) 



Atmospheric neutrino results give a Sm 2 from 10 2 to 
10 _3 eV 2 . Long base-line (L-B) reactor neutrino experi- 
ments |15j have been one of the most powerful ways to 
investigate I7 e — > v ^ neutrino oscillations (or, more gen- 
erally, v c -> 77 x oscillations). The CHOOZ [TM7| and 
PALO VERDE [TH] experiments utilized the high inten- 
sity and purity of the reactor core flux to achieve high 
sensitivity. 

The CHOOZ experiment had an average value of L /E ~- 
300 (L ~ 1km, E ~ 3MeV), an intense and nearly pure 
neutrino flavour composition (~ 100% V e ) and an inten- 
sity known to better than 2%. It could therefore make a 
definitive contribution to solving the problem of the at- 
mospheric neutrino anomaly. CHOOZ removed the possi- 
bility of explaining the atmospheric neutrino anomaly by 
v c Vfx oscillations and Super-Kamiokande showed that 
i/ft «-> v T caused the effect [T5) . 

The experiment was designed to detect reactor F e 's via 
the inverse /3-decay reaction 
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The signature is a delayed coincidence between the prompt 
c + signal (boosted by the two 511 — keV annihilation 7 
rays) and the signal from the neutron capture. The target 
material is a Hydrogen-rich (free protons) paraffin-based 
liquid scintillator loaded with Gadolinium, which was cho- 
sen due to its large neutron capture cross section and to 
the high 7-ray energy released after n-capture (<~ 8 MeV, 
well above the natural radioactivity). 

In this hnal paper we present a complete description of 
the experiment, calibration methods and stability checks, 
event reconstruction procedures and the Monte Carlo sim- 
ulation. 

Three different analyses, which have already been publis 
are more extensively discussed in this paper. The main 
one is based on all the available information: the mea- 
sured number of positron events as a function of energy, 
separately obtained from each reactor. It uses the two 
spectral shapes, as well as the absolute normalizations. 
The second result is based only on the comparison of the 
positron spectra from the two, different-distance reactors. 
This analysis is largely unaffected by the absolute value of 
the T> c flux, the cross section, the number of target protons 
and the detector efficiencies, and is therefore dominated by 
statistical errors. The sensitivity in this case is limited to 
Sm 2 > 2 • 1CT 3 eV 2 due to the small distance, AL = 116.7 
m, between the reactors. The explored (5m 2 , sin 2 26) pa- 
rameter space still matches well the region of the atmo- 
spheric neutrino anomaly. The third analysis is similar to 
the first, but does not include the knowledge of the abso- 
lute normalizations. 

Finally, some new remarks are presented, concerning 
the deduction of the confidence limits and the correct 
treatment of the systematic errors. 



2 The v c source and the expected interaction 
rate 

Nuclear reactors were the first sources used to search for 
neutrino oscillations and, in general, for systematic stud- 
ies of neutrino properties. In past experiments (i.e. up to 
the 1980s), the knowledge of the reactor neutrino flux and 
spectrum ( rj 10% accurate at that time) limited the sensi- 
tivity to oscillations. Oscillation tests were performed by 
comparing the neutrino event rate at different distances 
(up to ~ 100 m) from the reactor. This approach elimi- 
nated the systematic uncertainty associated with the ab- 
solute neutrino flux, often much greater than the experi- 
ment statistical error. Since then, our knowledge base of 
fission reactors, in particular of pressurized water reactors 
(PWR), has much improved, thanks to the previous reac- 
tor experiments and to direct studies of fission reactions 
on several elements. Present uncertainty on the neutrino 
flux and spectrum induces systematic errors on the event 
rate which are lower than 3%. It is therefore possible to 
perform experiments, at one fixed distance from a reactor 
only, relying on the adequate knowledge of the neutrino 
source. The reactor V e flux is perfectly isotropic and es- 
sentially with no contamination from other neutrino types, 



since the u c emission rate is lower than the one of v e by a 
factor > 10 5 and may therefore be discarded |20| . 



2.1 Description 

The CHOOZ nuclear power plant is located in the village 
of the same name, standing on the banks of the River 
Meuse in the north of France close to the border with 
Belgium. 

The power plant consists of two twin pressurized-water 
(PWR) reactors belonging to a recently developed gener- 
hetf'i^ * n France. The main innovation consists in an 
improved power yield (4.25 GWth, 1.45 GWe at full op- 
erating conditions), larger than any other PWR reactor. 
The first reactor reached full power in May 1997, the sec- 
ond in August 1997. Tab. summarizes data acquisition 
(from April 7, 1997 to July 20, 1998). The schedule was 



Table 1. Summary of the CHOOZ data acquisition cycle be- 
tween April 1997 and July 1998. 





Time (h) 


JWdt (GWh) 


Run 


8761.7 




Live time 


8209.3 




Dead time 


552.4 




Reactor 1 only ON 


2058.0 


8295 


Reactor 2 only ON 


1187.8 


4136 


Reactors 1 & 2 ON 


1543.1 


8841 


Reactors 1 & 2 OFF 


3420.4 





quite suitable for separating the individual reactor contri- 
butions and for determining the reactor-OFF background. 

The core of both reactors consists of an assembly of 
205 fuel elements bound to the socket plate of the reactor 
vessel. The vessel is filled with pressurized water (p = 155 
bars) at a temperature ranging from 280 °C at the en- 
trance to about 320 °C at the exit. The water, which acts 
as a neutron moderator and cooling element, circulates 
through four independent systems. Each of these involves 
a primary pump and a steam generator comprising 5416 
tubes immersed in the water of a secondary loop at a 
much lower pressure (56 bars) than primary loop pres- 
sure. As soon as the primary water passes through these 
tubes, the secondary water is vaporized and the steam 
produced travels to the turbine-alternator unit connected 
to the electric system [21"] . 

The water in the primary loop is slightly doped by 
Boron, a strong absorber of thermal neutrons. Boron ac- 
tivity yields information on the trend of the chain reaction 
in the core. 



2.2 Reactor power monitor 

Two methods were developed by the E.D.F. technical staff 
to monitor the total thermal power produced by each re- 
actor. The first one is based on the heat balance of the 
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steam generators. Measurements are made of parameters 
(water flow rate, vapour pressure and temperature in the 
secondary loop) needed to determine the amount of heat 
supplied to the steam generators. The resulting values are 
periodically available on a special computer network which 
records the data on an Excel file. The overall precision on 
the calculated thermal power is claimed to be 0.6%. 

The second set of thermal power data is provided by 
the external neutron flux measurement. This flux is ex- 
pected to be directly proportional to the fission rate in- 
side the core and, as a consequence, to the released ther- 
mal power. For each reactor, four different neutron de- 
tectors (one proportional counters plus three ionization 
chambers) are located at the opposite sides of the reac- 
tor vessel. This method is less precise than the other one 
because of the spread in the energy release per fission of 
the different fissile isotopes. Accuracy is estimated to be 
about 1.5%. However, this method has the advantage of 
operating continuously and is used simultaneously as a 
power monitor and sensor for the reactor safety system. 
Neutron detector outputs are fed to the computer network 
and written twice a day (or more frequently if the power 
ramps up or down) to another Excel file. The computer is 
also programmed to drive a direct current generator whose 
amplitude is proportional to the neutron detection rate. 
Such a signal was measured in the CHOOZ experiment by 
a CAM AC voltmeter across a 5J7 resistor (connected in se- 
ries to each loop). This provided a direct, instantaneous 
information about the thermal power for both reactors. 
The voltmeter calibration is shown in Fig. ^ where ther- 
mal power values, provided by E.D.F., are plotted versus 
the corresponding voltmeter data. The calibration param- 
eters obtained by a linear fit are also given. 




p0 = (-14743 ± 0.8) MW a 

pi = (14.837 + 0.0OSJ MW niV" f 



voltage (mV) voltage ImV] 

Fig. 1. Voltmeter calibration for reactor 1 (left) and 2 (right). 



are changed at the end of each cycle. The remainder is 
displaced towards the centre and new fuel elements are 
arranged in the outer part of the core, so as to get the fuel 
burning as uniformly as possible. For the start-up of the 
CHOOZ reactors, brand-new fuel rods were used. In or- 
der to reproduce the behaviour of reactors at equilibrium, 
where partially burned-up fuel is used, 137 less enriched 
elements are located at the centre (68 at 2.4% and 69 at 
1.8%, while the standard enrichment is 3.1%). A schematic 
map of the reactor core is drawn in Fig. 
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Region 1: 1.8 % enrichment 
Region 2: 2.4 % enrichment 

Region 3: 3.1 % enrichment 

Fig. 2. Schematic view of the fuel rods in the core for the first 
cycle of the CHOOZ reactors. The number of Boron poison 
rods assembled with each fuel element is also indicated. 



As shown in Fig. |3 the fuel assemblies of the most 
enriched elements are surrounded by a number of 12.7% 
Boron-doped steel rods, termed "poison" rods. The rods 
absorb the thermal neutron excess, thus accumulating 11 B, 
which has a negligible neutron absorption cross section, in- 
side the core. Therefore the neutron absorption power of 
the poison rods diminishes with time partly compensat- 
ing the fuel burn-up. The number of rods per fuel element 
varies according to the fuel loading as well as to the el- 
ement position in the core. This was taken into account 
when computing fuel evolution. 



2.3 Map of the reactor core 

The nuclear fuel consists of 110 T of uranium oxide tablets 
(diameter = 8.2 mm) enriched with 235 U and stacked in 
4 m long, 1 cm wide assemblies. The standard enrichment 
for this type of reactor is 3.1%. Each fuel element con- 
tains 264 assemblies. About 1/3 of the 205 fuel elements 



2.4 Fuel evolution 

The unit used to describe the aging of the fuel at nuclear 
reactors is the MWd/T, which measures the amount of 
energy per ton extracted from the nuclear fuel after its 
introduction into the reactor core. This quantity is called 
"burn-up" and is closely related to the fissile isotope com- 
position of the fuel. 
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For any PWR reactor, the procedure to compute the 
evolution of the fuel in the core needs daily information 
provided by the reactor technical staff. This includes: 



the daily cumulative burn-up, given as 



Mu J 



(3) 



W being the thermal power, t the time since the start 
of reactor operation and Mjj = 110.26 T the total 
amount of uranium in the core: 

the burn-up fti and the relative contribution to the 
power from the i-th fuel element, at several stages of 
the reactor cycle (an example is shown in Fig- EJ - 



P = 1000 MW d/T 



1 1.121 

, MOO , UfltL 



1.1V V. 
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Mflfl™- - 



i'u°§§ -' CTtTO^TO - Contribution to power 



» 9 9Q a l $j § 90 ^ Element burnup 

9(iri 7QQ 

Region index 

Fig. 3. Power distribution and burn-up values for the fuel 
elements in an octant of the CHOOZ reactor core at a certain 
step (/3 = 1000 of the first cycle). The contribution to the power 
of each element is normalized to have a mean value equal to 
one. 



These two inputs determine the daily burn-up of each fuel 
element. An example of the tables provided by E.D.F. for 
these quantities is given in Fig. E.D.F. also provides 
another set of tables (at several burn-up stages) in which 
the relative power contribution /? from the k-th fissile 
isotope for the i-th fuel clement is given. 

The number n\ of fissions per second for the i-th ele- 
ment for each isotope k can then be computed with the 
equation 



(4) 



where Ek is the energy release per fission for the fc-th iso- 
tope, whose values are listed in Tab. [21 Adding the contri- 
bution of all the fuel elements yields the average number 
Nk of fissions per second for the fc-th isotope. Contribu- 
tions from other fissioning isotopes, such as 236 U, 240 Pu, 



Table 2. Energy release per fission of the main fissile isotopes 
(from ref. 



isotope 


energy (McV) 




201.7 ±0.6 


238 tj 


205.0 ±0.9 


239p u 


210.0 ±0.9 


241 Pu 


212.4 ± 1.0 



242 Pu, etc. amount to less than 0.1% and are therefore 
neglected. 

To obtain the source spectrum, in addition to the av- 
erage fission rate Nk of each of the four isotopes, we need 
the corresponding differential neutrino yield per fission 
Sk(E u ). The next section explains how these spectra are 
evaluated. 



2.5 The expected neutrino spectra 

We used the so-called "conversion" approach which is the 
most reliable and recent method to determine the V e spec- 
trum at reactors. This method utilizes measurements of 
the electron spectrum emitted by a layer of fissile material 
activated by thermal neutrons. The experimental electron 
spectrum is then converted into the V e one |23U24U25| . The 
most recent and precise measurements were reported by 
Schreckenbach et al. |2fiU27U28| . In the latter case thin foils 
enriched with the main fissile nuclei (about 1 mg) were ex- 
posed to an intense thermal neutron flux (as 3 • 10 14 s _1 ). 
A high-resolution (Sp/p = 3.5 • 10~ 4 ) /3-spectrometer was 
used to measure the momentum of the emerging electrons. 

The (3~ spectrum for each fissile isotope is approxi- 
mated by the superposition of a set of hypothetical allowed 
branches with amplitude a% and end-point Eq\ 



Sp{Ep) 



iSa[Ep,EQ,Z(El) 



(5) 



where S% is the spectrum shape of the i-th branch and 
the summation is over the branches with end-point larger 
than Ep. The spectrum weakly depends also on the charge 
Z (average charge of the /3-decaying nuclei with end-point 
Eq) because of the Coulomb interaction in the final state. 
The measured electron spectrum is then used to determine 
the set of values {ai^E^} by means of a fit procedure. 
Eq. 10), with the introduction of the best-fit parameters, 
reproduces the measured spectrum to better than 1%. 

The f3~ spectrum for each individual hypothetical branch 
is then converted into the correlated v c spectrum under 
the assumption that both the electron and the antineu- 
trino share the total available energy E\. Thus, for each 
branch with end-point Eq, the probability of emitting an 
electron with energy Ep is equal to the probability of hav- 
ing a V e of energy Eq — Ep. Inserting the fit parameters 
into (JSJ, one obtains 



S v {Eu) = £ a t SU(Eo - E u ), El Z{E\ 



(G) 



2. The v c source and the expected interaction rate 



5 



These yields contain a normalization error of 1.9% stem- 
ming from the error on the neutron flux and from the 
absolute calibration uncertainty of the spectrometer. The 
conversion procedure also introduces a global shape un- 
certainty into the neutrino spectrum, beyond the inher- 
ent experimental errors. The main sources of this addi- 
tional uncertainty, ranging from 1.34% at 3MeV to 9.2% 
at 8 MeV, are the scattering in the nuclear charge distri- 
bution and the higher-order corrections (higher Coulomb 
terms and weak magnetism, for which an uncertainty of 
the order of the correction term itself was assumed) . 

This method was applied to obtain the neutrino yields 
of the 235 U, 239 Pu and 241 Pu fissions. The resulting spec- 
tra are presented in Fig. Unfortunately, no experimen- 
tal data is available for 238 U which cannot be fissioned by 
thermal neutrons. We must therefore rely on the theoret- 
ical predictions to estimate the contribution to the 
v e spectrum by the 238 U fission products. Although these 
predictions are less reliable than direct measurements, it 
should be noted that the contribution to the number of fis- 
sions, due to this isotope, is quite stable and never higher 
than 8%. Thus any possible discrepancy between the pre- 
dicted and the real spectrum should not lead to significant 
errors. 



of characteristics known at a few percent level. The sta- 
tistical accuracy obtained in these experiments makes it 
possible to use these results to discriminate between the 
existing models of reactor neutrino spectra. 

The Bugey 3 collaboration [32J measured the positron 
energy spectrum at 15 and 40 m from the reactor core 
and compared its data with the results of a Monte Carlo 
simulation using the neutrino spectrum models proposed 
by [27 ,33 34 . As can be seen in Fig. the data per- 
fectly fit with the measurements made at Institute Laue 
Langevin at Grenoble (ILL), whereas there is a lower com- 
patibility with other models. 
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Fig. 4. Neutrino yield per fission of the listed isotopes, 
determined by converting the measured f3 spectra |27I28| . 



2.6 Systematic uncertainties of the neutrino spectrum 

In past experiments at reactors, e.g. at Gosgen [HOj and 
Bugey more than 3x 10 4 neutrino events were recorded 
at several reactor-detector distances R. Since no evidence 
for oscillations was observed and rates were found to be 
consistent with a 1/i? 2 law, these experiments can be in- 
terpreted as a check that a reactor is a neutrino source 



Sehreckenbaeh et al. 



positron energy (MeV) 

Fig. 5. Comparison of Bugey 3 data with three different re- 
actor spectrum models. The error bars include only statistical 
uncertainties. The dashed lines are the quadratic sum of the 
quoted error of the models and the error due to the energy 
calibration. 



An improved measurement of the integral neutrino flux 
was performed in 1993. The measurement used an integral 
type detector previously employed at the Rovno nuclear 
power plant and subsequently moved to Bugey 36 . In 
that apparatus only neutrons were detected in 3 He coun- 
ters, while positrons were not. The apparatus was installed 
at 15 m from the reactor core. About 3 x 10 5 neutrino 
events were collected so that the reaction rate was deter- 
mined with 0.67% statistical accuracy. The neutrino event 
rate, n„, corresponds to a certain average fuel composition 
and is related to the cross section per fission Of and the 
number of target protons N p by 



1 Wi 



AttR 2 (E f ) 



(7) 



The fuel composition, the thermal power W, the average 
energy per fission (Ef) absorbed in the reactor core and 
the distance R were provided by the E.D.F. -Bugey tech- 
nical staff. The efficiency of the neutron detection e was 
carefully measured; the overall accuracy was estimated to 
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be 1.4%. The experimental result was then compared to 
the expected neutrino flux, which can be inferred by in- 
troducing into JJJ the neutrino spectra Sk obtained at 
ILL 26 , 27 28 , the cross section for the reaction J5J and 
the reactor parameters in Q. The expected cross section 
per fission for reaction is given by 

/>oo 

af p = / a{E v )S{E v )&E v 

r°a (8) 

= J2fk cr(E u )S k {E v ) dE u = V f k a k 
k J ° k 

where f k refers to the contribution of the main fissile 
nuclei to the total number of fission, S k to their corre- 
sponding v c spectrum and a k to their cross section per fis- 
sion. The result of the measurement was <r™ eas = 5.752 x 
10 _19 barns/fission± 1.4%, which perfectly agrees with the 
expected value (af p = 5.824x 10- 19 barns/fission±2.7%), 
and is twice as accurate as the predictions based on the 
knowledge of the neutrino spectrum. 

We could therefore adopt the conversion procedure for 
the shape of the neutrino spectra but normalize the to- 
tal cross section per fission to the Bugey measurement, 
i.e. , after taking all the different reactors conditions into 
account. 

In figurelHJleft) the ILL cross section JHJ and the com- 
bined (ILL+Bugey) cross section obtained for the CHOOZ 
reactors are plotted vs. the reactor burn-up. The average 
ratio of the two curves amounts to 0.987. By combining 
the uncertainty on the neutrino spectra, on the cross sec- 
tion for the reaction J5J) and on the fission contributions 
fk (which are of the order of 5%), we obtained the rel- 
ative error on the neutrino detection rate, as a function 
of the fuel burn- up. As shown in Fig. fright), the aver- 
age error decreases from 2.4% (ILL data alone) to 1.6% 
(ILL+Bugey). Other minor sources of errors come from 
the residual neutrino emission from long-lived fission frag- 
ments (dealt with in the next Section). However the reac- 
tor source can be considered to be known at a 2% level. 



2.7 Neutrino spectrum time relaxation and residual 
neutrino emission 

Another possible source of uncertainty of the neutrino flux 
is related to the residual emission due to the (3~ decay 
of long-lived fission fragments. Taking this further contri- 
bution into account, the linear relation J7J between the 
prompt V c interaction rate and the current thermal power 
W no longer holds. Nevertheless, the V e spectra deter- 
mined at ILL were derived after about 1.5 d exposure time, 
so that neutrinos from fission fragment decays of longer 
life are not included. The expected neutrino rate based 
on this model may thus be underestimated with respect 
to the experimental data. Fortunately the maximum neu- 
trino energy is above the reaction (J2J threshold only in a 
few cases since, as can be expected, the longer the lifetime, 
the lower the Q-value of the decay 1 . This effect has been 
evaluated by using the cumulative yields of the known 
long-lived fission fragments |37]; the results for 235 U and 
239 Pu are summarized in Tab.|21 In particular, the reaction 



Table 3. Time evolution of neutrino spectra relative to infinite 
irradiation time (from |37)h 



E v (MeV) 




!Kfl)p u 


1.5 


10 4 s 


1.5 d 


10 y 


10 4 s 


1.5 d 


10 y 


0.837 


0.946 


0.988 


0.861 


0.949 


0.990 


2 


0.897 


0.976 


0.992 


0.904 


0.968 


0.986 


2.5 


0.925 


0.981 


0.990 


0.939 


0.975 


0.986 


3 


0.963 


0.997 


1.000 


0.967 


0.989 


0.993 


3.5 


0.967 


1.000 


1.000 


0.979 


0.997 


1.000 



cross section 07 computed by using © is probably lower 
than the effective one by « 0.3%. This systematic shift 
affects the accuracy on the reaction cross section. We will 
then assume an overall and conservative 1.9% uncertainty 
on the integral neutrino rate. 




P(MW d/ton) P(MW d/ton) 

Fig. 6. Comparison of the combined (ILL+Bugey) reaction 
cross section with the ILL cross section (left) and their relative 
error (right) as a function of the first CHOOZ reactor cycle 
burn-up. 



2.8 The inverse beta-decay reaction 

The detection of reactor antineutrinos is usually based on 
the reaction J5J. This is the most suitable process since: 

- it has the highest cross section in the energy range of 
reactor antineutrinos: 

- provides a convenient time correlated pair of positron 
and neutron signals, which allows us to reject most of 
the background. 

The antineutrino and the positron energy are related by 

£W C - E c+ + (M n - A/ p ) + 0(E v jM n ), (9) 

1 This effect might be much more relevant in the case of 
experiments looking for neutrino elastic scattering interactions 
(i.e. , measurement of the neutrino magnetic moment) and 
needs a more careful treatment. 
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where the infinitesimal term corresponds to the neutron 
recoil. Thus a measurement of the positron energy al- 
lows an accurate determination of the energy of the in- 
coming antineutrino. The threshold for the reaction J5J 
is 1.804 MeV, equal to the nucleon mass difference plus 
the positron mass. In the low energy limit, the cross sec- 
tion for the reaction (J2J) may be written as a function of 
outgoing positron energy as follows: 

2n 2 h 3 . , . 

a(E e +) = =- p c +E c +(l + Orad + 5wm) (10) 

The transition matrix element has been expressed in terms 
of the free neutron decay phase-space factor / = 1.71465(15) 
|58") and lifetime r„ = (886. 7± 1.9) s Two higher-order 
correction terms (both of 1% order of magnitude) are also 
included: 

(i) a radiative correction of the order of a, including an 
internal bremsstrahlung contribution, which can be ap- 
proximated by 



Srad(E e+ ) = 11.7 x 10-%E e+ - m e ) 



\-0.3 



(11) 



with the positron energy expressed in MeV 0U] an d 
(ii) a correction for weak magnetism, arising from the dif- 
ference [i — [i n — Hp = — 4.705890(2)/^ between the 
anomalous magnetic moment of the neutron and the 
proton 



(E e+ ) = -2 



//A 



1 + 3A 2 



(E e+ +A(3p e+ )/M p , (12) 



where A = gA/gv — 1.2601 ± 0.0025 is the ratio of 
axial-vector and vector coupling constants and A is 
the nucleon mass difference |41| . 

The knowledge of the cross section is therefore much more 
accurate than the one of the V c spectrum, the major lim- 
itation being related to the uncertainty on r„ (whose rel- 
ative error is at most 0.3%). 



2.9 Simulation of the v c spectrum 

The expected neutrino spectrum was obtained by means 
of a Monte Carlo simulation of the reactor core in which 
we folded all the ingredients described in the previous sec- 
tions, that is: 

— the daily variation of the flux due to the fissile elements 
burn-up and to the live time of the apparatus; 

— the contribution of the individual fuel elements accord- 
ing to their position inside the reactor core; 

— the individual contributions of the different fissile ele- 
ments in each fuel element to the neutrino flux. 

Fig.[7|shows the neutrino spectrum obtained at CHOOZ 
on the start-up day (J3 = 0) and at an intermediate step 
(j3 = 7000) of the first reactor cycle. Due to the decrease 
of the 235 U concentration, a reduction of the neutrino 
interaction rate is observed and a softening of the neu- 
trino spectrum is expected. The energy spectrum of the 
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Fig. 7. Comparison between the neutrino spectrum at the 
beginning and during the first cycle of the CHOOZ reactors. 



positrons coming from is essentially the antineutrino 
spectrum shifted in energy and weighted by the cross sec- 
tion l|10|) . So, following |0, each positron is assigned a 
weight given by 



S e +(T e -i 



N, 



And 2 



?-a{E v )S v {E v ), 



(13) 



where d is the distance from the neutrino production point 
to the detector and the positron kinetic energy T e + is given 
by {HI. Fig. |H1 shows the positron yield obtained by gener- 
ating the neutrino spectra drawn in Fig. in one day of 
data taking with both reactors at full power. Although the 
generated neutrino luminosity is the same, the decrease of 
the positron yield with the reactor operating time is evi- 
dent. The evolution of the positron spectrum must be fol- 
lowed very accurately in order to reproduce the hardware 
threshold effects on the positron detection. As an imme- 
diate consequence, also the integral neutrino interaction 
rate is expected to vary significantly during the reactor 
fuel cycle. A decrease of about 10% has been forecast for 
the cross section per fission (which is linear with the in- 
teraction rate, according to 0) during the first cycle of 
the CHOOZ reactors, as shown in Fig. [5] The measured 
neutrino rate as a function of the burn-up will be shown 
and compared to the expected behaviour, under the no- 
oscillation hypothesis. 



3 The experiment 
3.1 The site 

The detector was located in an underground laboratory 
about 1km from the neutrino source (see Fig. llOfl . The 
300 MWE rock overburden reduced the external cosmic 
ray muon flux by a factor of ~ 300 to a value of 0.4 m~ 2 s -1 , 
significantly decreasing the most dangerous background, 
which is caused by fast neutrons produced by muon-induccd 
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Fig. 8. Positron spectra at the startup and during the first 
cycle of the CHOOZ reactors at maximum daily neutrino lu- 
minosity. 
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Fig. 9. Cross section per fission as a function of the reactor 
burn-up. The contribution of each fissile isotope is also shown. 



nuclear spallations in the materials surrounding the detec- 
tor. This cosmic ray shielding was an important feature of 
the CHOOZ site. As shown in Fig. ^] the rock shielding 
preserved the signal to noise ratio of previous reactor ex- 
periments, in spite of a reduction by a factor ~ 100 of the 
neutrino flux due to the larger distance from the reactors. 
The characteristics of the site were thoroughly assessed 
before and after the set-up of the experiment. We mea- 
sured the cosmic ray flux and angular distribution and the 
results were compared with predictions based on the rock 
structure above and around the site. The natural radioac- 
tivity background, the radon level in the tunnel and the 
intensity and orientation of the local terrestrial magnetic 
field were also measured. A geodesic survey was performed 
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Nuclear Power Station 
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Chooz Underground Neutrino Laboratory 
Ardennes, France 



Fig. 10. Overview of the experiment site with indication of 
the source-detector distance and rock overburden. 



to determine the exact location of the underground exper- 
iment with respect to the two reactors and its orientation. 
All the information obtained in the preliminary studies 
were used to guide improvements in the detector and in 
the detector simulation. 

The cosmic ray measurements were made for several 
days with a system of six Resistive Plate Chambers (RPC), 
each of area lxl m 2 . The comparison of the experimental 
and expected angular distributions was fairly good, but 
some discrepancies persisted. A further geological study 
revealed the existence of several very high density rock lay- 
ers (3.1 g/cm 3 the normal density being 2.8 g/cm 3 ), whose 
positions and orientations fully explained the observed ef- 
fects. The natural radioactivity spectrum was measured, 
at the position where the experiment had to be installed, 
by a 3 x 3 inches 2 Nal crystal. Natural radioactivity was 
rather high and the spectrum shows the normal natural ra- 
dioactivity lines but also some artificial radioactivity con- 
tributions (see Fig. 1X1)1 . 

The magnetic field was measured by a rotating coil, 
pressurized-air device, and was found to be B\\ = 0.178 ± 
0.007 G and B± = 0.388 ± 0.007 G. 

3.2 The detector 

The detector (Fig.JXjJ was installed in a welded cylindrical 
steel vessel 5.5 m in diameter and 5.5 m deep. The inter- 
nal walls of the vessel were painted with high-reflectivity 
white paint. The vessel was placed in a pit 7 m in diame- 
ter and 7 m deep. To protect the detector from the natural 
radioactivity of the rock, the steel vessel was surrounded 
by 75 cm of low radioactivity sand (Comblanchien from 



3. The experiment 



9 



RADIOACTIVITY SPECTRUM AT CHOOZ SITE 



7 



10 



A ! " Bi B7 

1 \ : 13, Cs 




I I : 


A, I 2,4 


Si 




~ 2IJ Pb Ssp 




208™ 1 
208™ 11 = 


M Co 


1 j 2 " Bi 






214 Bi 


_ 232 Th Family \ 

238 U Family \ 


Artificial Radioactivity \ 

i I...." ill 



ENERGY (KcV) 

Fig. 11. Natural radioactivity spectrum recorded by a Nal 
crystal in the well hosting the detector. 
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Fig. 12. Cosmic muon flux compared to the neutrino flux at 
the different underground experimental sites. In the CHOOZ 
case the lower neutrino flux is compensated by the reduction 
of the muon flux. 



Burgundy in France) and covered by 14 cm of cast iron. 
The detector comprised three concentric regions: 

— a central 5-ton target in a transparent Plexiglas con- 
tainer (total mass = 117 kg) filled with a 0.09% Gd- 
loaded scintillator ( "Region I" ) ; 

— an intermediate 17-ton region (70 cm thick) equipped 
with 192 eight-inch PMT's (15% surface coverage, ~ 
130 photoelectrons/MeV), used to protect the target 
from PMT radioactivity and to contain the gamma 
rays from neutron capture ("Region II"); 

— an outer 90-ton optically separated active cosmic-ray 
muon veto shield (80 cm thick) equipped with two rings 
of 24 eight-inch PMT's ("Region III"). 

The apparatus was conceived as a liquid scintillator low 
energy, high-resolution calorimeter. The detector geome- 
try (a central volume of scintillator surrounded by photo- 
multipliers) is common to the Borexino, LSND and SNO 
detectors. The detector was simple and easily calibrated, 



while its behaviour could be well checked. Six laser flashers 
were installed in the three regions together with calibra- 
tion pipes to allow the introduction of radioactive sources. 
The detector could be reliably simulated by the Monte 
Carlo method. 

The target region contained a Gd-loaded liquid scin- 
tillator. The neutrino detection was based on the delayed 
coincidence between the prompt positron signal generated 
by reaction (0), boosted by the annihilation 7-rays, and 
the signal associated with the 7-ray emission following the 
neutron capture reaction 

n + Gd -> Gd* -> Gd + ^ 7l (14) 

i 

The choice of a Gd-doping was to maximize the neutron 
capture efficiency; Gadolinium has the highest thermal 
neutron cross section. Moreover, the large total 7-ray en- 
ergy (» 8 MeV, as shown in Tab. 0} easily discriminates 
the neutron capture from the natural radioactivity, whose 
energy does not exceed 3.5 MeV. 

Table 4. Abundances and thermal neutron capture cross sec- 
tions for the Gd isotopes. 



Gd 




Abundance 


Cross section 


Relative 


isotope 


(KeV) 


(%) 


(barns) 


intensity 


152 


6247 


0.20 


735 


3 ■ 10~ 5 


154 


6438 


2.18 


85 


3.8- 10~ b 


155 


8536 


14.80 


60900 


0.1848 


156 


6360 


20.47 


1.50 


6 • 10~ b 


157 


7937 


15.65 


254000 


0.8151 


158 


5942 


24.84 


2.20 


1.1 ■ 10~ b 


160 


5635 


21.86 


0.77 


3 • 10~ b 



Region II was filled with an undoped high-flash point 
liquid scintillator. It provided a high-efficiency contain- 
ment of the e.m. energy deposit; this was higher than 99% 
for positrons from T7 c -interactions in Region I. The con- 
tainment of the 7-rays due to the neutron capture on Gd 
was (on average) slightly lower than 95% for an energy de- 
posit E > 6 MeV. The intermediate volume was bounded 
by the "geode" , an opaque plastic structure serving as a 
support for the 192 inward-looking photomultiplier tubes 
(PMT from now on). 

The outer volume, also filled with the undoped scin- 
tillator of Region II, was the "Veto" region (Region III). 
An additional 48 PMT's, arranged in two circular rings 
located at the top and the bottom of the main tank, de- 
tected the scintillation light associated with through-going 
cosmic muons. The Veto signal was used to tag and reject 
this major background source. The outer scintillator layer 
was also thick enough to shield the neutrino target against 
the natural radioactivity from the surrounding materials. 

The inner detector volume was separated from Region 
II by a transparent 8 mm-thick vessel, a vertical cylindrical 
surface closed by two hemispherical end-caps. The outer 
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radius of the cylinder and of the end-caps was 90 cm, the 
height of the cylinder was 100 cm; the inner volume was 
5.555m 3 , while the mass was 150kg (empty). The vessel 
was made of an acrylic polymer (Altuglass), chosen for 
its excellent optical and mechanical properties and for its 
chemical resistance to aromatic compounds in the scin- 
tillator. The upper part of the vessel was fitted with a 
chimney (diameter = 70 mm) to allow passage of fill- 
ing pipes, calibration sources, temperature and pressure 
sensors. 

The geode had the same shape of the target vessel, but 
a larger size; the cylinder height was the same, whereas the 
outer radius was 160 cm. The volume between the geode 
and the target was 19.6 m 3 . The geode surface (a drawing 
of which is shown in Fig. I13[) had a total area of 42 m 2 
segmented into 32 panels; each panel was equipped with 
6 8" PMT's detecting the scintillation light produced in 
Regions I and II. The global PMT coverage was then 15%. 
Unlike the acrylic inner vessel, the geode was opaque so as 




Fig. 13. Mechanical drawing of the detector; the visible holes 
on the geode are for the PMT housing (from \4'2]). 



to optically shield the inner regions from the veto scintilla- 
tion light. The external surface was white-coated in order 
to enhance the light collection in Region III and improve 
the Veto rejection efficiency, while the black inner surface 
reduced light reflections, which could degrade the vertex 
resolution in the detector |43|. 



3.3 The liquid scintillators 

About 5 tons of Gd- loaded scintillator and 107 tons of 
unloaded scintillator were used in the experiment. The 
main properties of the two scintillators are listed in Tab. [5] 



Table 5. Main properties of the liquid scintillators used in the 
experiment. 





Gd-loaded 


unloaded 


Chemical content: 






basic 


Norpar-15 


Mineral oil 




(50% vol.) 


(92.8% vol.) 


aromatics, alcohols 


IPB+hexanol 


IPB 




(50% vol.) 


(7.2% vol.) 


wavelength shifters 


p-PTP+bis-MSB 


PPO + DPA 




(1 g/1) 


(1.5 g/1) 


Atomic mass 






composition: 






H 


12.2% 


13.3% 


C 


84.4% 


85.5% 


Gd 


0.1% 




others 


3.3% 


1.2% 


compatibility 


acrylic, Teflon 


density (20 °C) 


0.846 g/ml 


0.854 g/ml 


Flash point 


69 °C 


110°C 


Scintillation yield 


5300 ph/MeV (35% of anthracene) 


Optical attenuation 


4 m 


10 m 


length 






Refractive index 


1.472 


1.476 


Neutron capture 


30.5 /is 


180 /xs 


time 






Neutron capture 


~ 6 cm 


~ 40 cm 


path length 






Capture fraction 


84.1% 




on Gd 







The solution of the Gadolinium salt Gd(NOs)3 in hex- 
anol as well as the mixing of the basic and aromatic com- 
pounds was performed in a dedicated hall close to the en- 
trance of the underground tunnel. The amount of Gadolin- 
ium (0.09 % in weight) was chosen to optimize neutron 
capture time and neutron detection efficiency. The mea- 
sured values for the average capture time, path length and 
capture efficiency are listed in Tab.JSJ A higher concentra- 
tion would have required more alcohol, which could have 
lowered the high flash point of the solution below the limit 
imposed by safety regulations. Moreover, as we shall see, 
the presence of the nitrate ions in solution progressively 
deteriorated the optical properties of the scintillator. A 
higher concentration would have further compromised the 
chemical stability of the scintillator. 

A fundamental quantity for normalizing the neutrino 
event rate is the number of free protons in the target. 
An accurate evaluation of this number relies on precise 
measurements of the density and of the Hydrogen con- 
tent of the Gd-loaded scintillator. The Hydrogen content 
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was determined by combustion of scintillator samples; the 
presence of volatile elements (chiefly various types of al- 
cohol) made this measurement particularly difficult. The 
value listed in Tab. El derives from averaging the dozen 
measurements performed on scintillator samples at differ- 
ent laboratories | 44U45 | . The overall relative uncertainty 
on the number of protons is 0.8%. 

The scintillator transparency showed a degradation over 
time, resulting in a slight decrease of the photoelectron 
yield. The most probable cause is thought to be the ox- 
idation by the nitrate ion. On one side, this forced us 
to repeatedly check the attenuation length in the detec- 
tor by using calibration sources all along the acquisition 
period. The procedure followed is described in j EH In 
parallel, we made a systematic laboratory study of the 
chemical and optical stability of liquid Gd-loaded scin- 
tillators, which is also of interest for future low-energy 
neutrino detectors. The main optical properties were de- 
termined with both spectrophotometric and fluorimetric 
measurements. A Varian Cary 2200, UV-vis, 10 cm cell 
double-beam spectrophotometer was used to measure the 
absorbance of scintillators as a function of the beam light 
wavelength in the 350 600 nm range, the absorbance 
being defined as 



log! 



Ia 



(15) 



where in,/ are the intensity respectively of the incident 
and the emerging beams. Several paired cells, with Suprasil 
quartz windows were used for these measurements. The 
measurement of attenuation lengths of a few meters by 
means of 10 cm cells required great care in the stability 
control of the apparatus as well as an accurate cleaning of 
the cell optical surfaces. Corrections were applied to take 
into account the multiple reflections inside the cell due to 
the different refraction indices of air, quartz and scintilla- 
tor; the overall systematic uncertainty on the absorbance 
(including the spectrometer stability and the planarity of 
the quartz windows) was proved to be 0.5%. Examples of 
these measurements obtained with a laboratory-blended 
scintillator sample (using the same composition as the Re- 
gion I scintillator) are shown in Fig. 1141 (left). 

An overall attenuation for the scintillation light was 
obtained by folding the measured attenuation at different 
wavelengths with the fluorescence spectra of the phospho- 
res and the PMT photocathode sensitivity. The light at- 
tenuation as a function of the path is shown in Fig. 1141 
(right). Apart from the first ~ 20cm, needed for the ab- 
sorption by the second shifter (bis-MSB), attenuation is 
best described by an exponential decrease giving an "effec- 
tive" length. Measurements performed on the laboratory 
samples at different aging stages reproduced the atten- 
uation length values obtained for the target scintillator 
within 15%. 

The time variation of the light attenuation length in 
the detector is best reproduced by the function 



d5T 



Fig. 14. Attenuation length vs. wavelength for the Gd-loaded 
scintillator (of the same type of that used at CHOOZ) at dif- 
ferent aging stages (left) and scintillation light attenuation vs. 
path (right). 



which accounts for the observed exponential decrease of 
the number of photoelectrons with time. Here v is propor- 
tional to the velocity of the chemical reactions responsible 
for the scintillator deterioration. The reaction kinetics is 
known to depend on temperature according to an expo- 
nential law: 



w (T) = Uo /(r) = W (T - To)/loocl 



(17) 



\(t) = 



An 



1 + vt 



(16) 



where the index labels the values at room temperature. 
It is also known that reactions occurring at time scales 
~ 1 s are accelerated by a factor ~ 2 if the temperature is 
raised by 10 °C from room temperature. 

We attempted to accelerate the aging effects by heat- 
ing different samples at 60, 70 and 80 °C. As shown in 
Fig. we found that an increase of 10 degrees in the 
temperature corresponded to an acceleration factor a ~ 3 
(instead of 2) in the aging rate. This discrepancy can be 
explained by simple thermodynamical arguments if (as in 
this case) such a reaction develops at room temperature 
on a time scale of a few months. The estimated value (at 
T = 20 °C) was v = (3.8 ± 1.4) • lO^d" 1 , which was 
acceptable for a 1-year data-taking duration experiment 
like CHOOZ. This prediction turned out to be close to the 
value ((4.2 ± 0.4) • 10 _3 d _1 ) obtained by direct measure- 
ments in the detector. 



3.4 The photomultiplier system 

The PMTs provided both the charge and the time infor- 
mation needed to determine the event position and the 
energy deposit. The average photoelectron (pe) yield as- 
sociated with a typical positron signal (E w 3MeV) was 
« 300, which corresponded to an average « 1.5 pe at each 
PMT. The number of hit PMTs was w 80% of the total. 
As the PMTs had to work in a single pe regime, the most 
important feature for tube selection was the single pe res- 
olution, i.e. , the peak to valley ratio. Other important 
requirements were: low radioactivity and dark noise, high 
gain, large photocathode surface, low sensitivity to mag- 
netic fields, good (but not ultra-fast) time properties. The 
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Fig. 15. Acceleration of the scintillator aging rate as a function 
of the temperature. 



(iii) the relative light sensitivity (proportional to the pho- 
tocathode quantum efficiency x pe collection efficiency); 

(iv) the single pe pulse height spectrum and its peak to 
valley (P/V) ratio; 

(v) the PMT time characteristics (time jitter and "walk" 
effects). 

The results were used to define the optimal working con- 
ditions and decide a proper geometrical arrangement of all 
the PMT's in the detector. Since the light level of neutrino 
events in CHOOZ corresponded to a few photoelectrons 
at each PMT, it is at this regime that the listed param- 
eters were determined. The best suited source was the 
Hamamatsu light pulser PLP-02 with a laser diode head 
SLD-041 emitting at 410 nm a wavelength close to the 
maximum photocathode sensitivity and to the maximum 
bis-MSB emission wavelength. Details of the features of 
this source, the coupling to the optical bench and to PMT 
housing structure can be found in @7J. The determination 
of the PMT relative sensitivity was cross-checked by us- 
ing two additional "passive" light sources: a disk of NE110 
scintillator activated by a low intensity 60 Co and a 241 Am 
a-source coupled to a Nal crystal. All our sensitivity mea- 
surements were found to be consistent with each other and 
with the quantum efficiency as measured by EMI. 



EMI 9536KA B53 |U turned out to be the best PMT 
for our needs. The high amplification (« 10 , correspond- 
ing to w 30 mV for a 1-pe pulse) matched the character- 
istics of the front-end electronics; the transit-time jitter 
(« 8 ns FWHM at 1-pe level) was comparable with the 
scintillator decay time (rj 7 ns). The PMT envelope was 
made of an ultra low-radioactivity glass (B53) at our re- 
quest. PMTs were power supplied by the CAEN SY527 
HV system (positive anode, grounded cathode), equipped 
with 16-channel cards A734P, providing a maximum of 
3 kV; the HV channels were CAMAC controlled via the 
CAENET card C117B. The high voltage was distributed 
to the dynodes through a "tapered bleeder" progressive 
divider, with interdynode voltages progressively increas- 
ing in the last dynode stages. This divider was preferred 
to the linear one (with equal interdynode voltages), since 
it was proved to provide enhanced linearity at a still ac- 
ceptable high voltage for the required gain. The divider 
consisted in a circular printed circuit board soldered to 
the PMT socket pins and glued to the back of the socket 
itself. The divider and the socket pins were enclosed in a 
Plexiglass cylinder to prevent contact of liquid scintillator 
with all the components. 

A PMT test facility 07J was designed to determine 
the relevant parameters of all the PMTs. The measure- 
ment program allowed the quantities listed below to be 
determined and thus to compare the results with the ones 
contained in the PMT test card given by the producer: 

(i) the operating voltage corresponding to 30 mV pulse 
(on 50 Q) for a single pe; 

(ii) the PMT noise level; 



3.5 Detector simulation 

The Monte Carlo (MC) simulation of the detector is based 
on a Fortran code linked with the CERN GEANT 3.21 
package The program allowed us to simulate the de- 
tector response for different particles, in particular for the 
positron and the neutron following an V c interaction. The 
use of GEANT routines had a two-fold objective: 

- define the detector geometry and the physical param- 

eters of the different materials (scintillators, acrylic, 
steel) in the detector and 

- track positrons, electrons and photons in the detector 

and evaluate the energy deposit in the scintillator. 

A specially written routine was used to track secondary 
particles below 10 keV, the minimum cut-off energy al- 
lowed by GEANT, down to 1 keV, in order to better ac- 
count for the scintillation saturation effects arising at low 
energies. The routine calculates the range corresponding 
to the particle energy and divides this range into 10 steps; 
the ionization loss is evaluated at each step, and the parti- 
cle is tracked until its residual energy reaches 1 keV. Each 
energy deposit is then converted into the scintillation pho- 
ton yield; ionization quenching effects are taken into ac- 
count according to Birk's law, with j3 = 1.3T0 -2 cm/ MeV 
derived from the measurements performed on the similar 
MACRO liquid scintillator g2J. 

Scintillation photons are isotropically generated and 
tracked up to the PMT surface. Effects due to light atten- 
uation by the scintillators, due to reflection on the PMT 
glass, to the light refraction at the target boundary and 
to absorption on opaque surfaces (such as the Geode and 
the calibration pipes) are also simulated. The probability 
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of detecting a photon hitting the PMT surface is weighted 
according to the quantum efficiency. 

The MC code also includes the charge response of the 
digitizing electronics; the charge output for each event is 
formatted as for real data, which allowed us to reconstruct 
MC and real events by means of the same algorithm. The 
evaluation of the neutrino detection efficiencies strongly 
relies on MC evaluations; it was therefore crucial that the 
MC output was as close as possible to the data. The most 
important difference between the real and simulated ge- 
ometry concerns the shape of the PMT's, assumed to be 
flat and lying on the geode surface (whereas the PMT 
glass is hemispherical and protrudes towards the centre 
by a few cm). As we will see, this approximation (needed 
to avoid a large waste of computing time) is also responsi- 
ble for the fragility of the reconstruction algorithm in the 
case of events close to the PMTs. 

User-defined routines j50| were introduced in the main 
MC code to follow the neutron moderation and capture 
processes. Tabulated cross section values were used for 
both elastic scattering and capture on Gd, H, C, O, Fe |51| . 
At each step (whose length depends on the total cross sec- 
tion at the neutron velocity) a decision is taken if either 
scattering or capture processes (or neither) arise. In the 
case of captures, 7-rays are emitted in cascade according 
to the de-excitation scheme of the nucleus involved; each 
photon is then tracked by GEANT routines. Captures on 
Hydrogen are followed by the emission of one 2.2 MeV 
7-ray. The case of Gadolinium is harder to manage; the 
excited nucleus can decay to the ground state through 
a series of intermediate excited levels, both discrete and 
continuously distributed. Each capture event typically re- 
leases about three 7-rays with a total energy depending 
on which Gd isotope captures the neutron. The informa- 
tion used about the energy levels involved and the relative 
decay branches are taken from |52| . 



4.1 The first-level trigger 



The electronic sum of the signals corresponding to photo- 
electrons detected by all PMT's gave a first rough mea- 
surement of the energy released by any event. A Monte 
Carlo simulation of electrons uniformly generated in the 
detector (whose results are summarized in Fig. ll6|) showed 
that the collected charge is linearly dependent on the en- 
ergy deposit and almost independent of the position within 
Region I. In Region II this changes and the total number 
of photoelectrons for events within 30 cm from the PMT's 
can be as much as 10 times larger than for events at the 
centre. So, a trigger exclusively based on the charge infor- 
mation would not have rejected lower energy events gen- 
erated close to the PMT surface. 
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Fig. 16. Number of detected photoelectrons QSUM (left) and 
number of hit PMT's NSUM (right) for 1 MeV electron events, 
as a function of the distance from the detector centre. 



4 The trigger 

The electronic system for the CHOOZ experiment col- 
lected and processed PMT output signals to select events 
due to neutrino interactions and to separate them from 
the background. The charge and time information pro- 
vided by the PMT's were used to reconstruct the energy 
and position of each event. 

The signature of a neutrino interaction relied on the 
delayed coincidence of two signals within a 100 fj,s time 
window: the former, due to the positron, has an aver- 
age energy of about 3 MeV; the latter, due to the de- 
layed neutron capture on Gadolinium, has a total 7-ray 
energy w 8 MeV. The identification of this event com- 
bination was therefore based on a two-level trigger sys- 
tem; the first level selected events with an energy release 
> 1.3 MeV (Lllo) or > 3 MeV (Llhi); the occurrence of 
the sequence (Lllo, Llhi), within the 100 [is time window, 
generated the second level which triggered the data digi- 
tization and recording by the on-line system. 



A "topological" cut was then applied by also requir- 
ing a minimum number of hit PMT's. Fig. llGl clearlv shows 
that this extra condition preferentially rejected events close 
to the PMT's. 

The first-level (LI) trigger was generated when both 
the charge and the PMT multiplicity fulfilled the above 
criteria. A simple schematic of the trigger circuit is pre- 
sented in Fig. El The signals from each PMT were fed into 
the front-end electronics made of fan-in/out modules de- 
veloped for the purpose. These modules provided a linear 
sum of the input signals (QSUM) whose amplitude was 
proportional to the total number of photoelectrons. The 
PMT multiplicity signal (NSUM) was obtained in a sim- 
ilar way. A copy of the PMT signals was fed into Lecroy 
LRS 4413 discriminators which had their threshold ad- 
justed to 15 mV (« one half of the single photoelectron 
amplitude). The linear sum of these discriminator out- 
puts, whose amplitude was proportional to the number 
of channels over threshold on each discriminator board, 
formed the NSUM signal. The LI trigger was finally as- 
serted when both the resulting QSUM and NSUM signals 
exceeded the preset thresholds. 
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Fig. 17. First-level trigger scheme. Both the number of photo- 
electrons (QSUM) and the number of hit PMT's (NSUM) are 
required to fulfil a certain threshold condition. 

4.2 The second level trigger 

The second-level (L2) condition (the one identifying the 
neutrino interactions) triggered on the occurrence of two 
events satisfying the LI condition within a 100 fis wide 
time window. This width represented a trade-off between 
a high neutrino detection efficiency and a high rejection 
of the accidental background. The positron-neutron delay 
followed an exponential distribution with an average r ~ 
30 /is, therefore only 5% of the neutrino events were missed 
by the delayed coincidence window. 

The L2 trigger logic is somewhat more complicated 
than the simple scheme just described. There were two 
different conditions for the fulfilment of the LI trigger, 
corresponding to two different energy thresholds: the low 
condition (Lllo) roughly corresponded to 1.3 MeV and 
had a rate w 130 s _1 , the high condition (Llhi) corre- 
sponded to w 3 MeV and had a rate ~ 30 s -1 . The higher 
threshold rejected the large background due to the natural 
7-radioactivity (which practically ends with the 2.6 MeV 
208 T1 line) and was sufficiently low to detect signals due 
to the neutron capture on Gd with almost full efficiency. 

The L2 trigger was vetoed if some activity was present 
in Region III; this selection rejected events associated with 
cosmic rays. The VETO condition was satisfied whenever 
VSUM, the sum of the signals of the 48 PMT's located 
in Region III, was over a preset threshold. If no VETO 
condition occurred during a 1 ms time interval preceeding 
L2, the L2 triggered: the acquisition of further secondary 
particles went on for an extra-time of 100 /is and then 
stopped. The on-line processor took « 80 ms to read all 
the electronics and to record the data on disk. During 
this time the computer busy condition was on and the 
acquisition of any further event was disabled. 

Although a neutrino interaction corresponded to an 
Lllo-Llhi or Llhi-Llhi time sequence, the trigger logic en- 
abled the acquisition of events with an Llhi-Lllo sequence 
as well, allowing a systematic study of the accidental back- 
ground. 

Finally, although the acquisition rate (« 0.15 s _1 ) was 
rather low, it was large compared to the expected neutrino 
rate at full reactor power (« 30 events d _1 ). The typical 
event size was roughly 30 Kbytes; the acquisition system 



could handle a daily amount of data in the order of 0.5 
Gbytes. This capability was achieved through fast readout 
electronics and large storage devices. 



4.3 Trigger for calibration runs 

Acquisition schemes other than the normal neutrino trig- 
ger were required for calibrations or particular test runs. 
The trigger logic could be changed at the beginning of 
each run by loading an appropriate trigger table, a combi- 
nation of the level trigger, veto and busy signals assembled 
by a dedicated logic unit Caen C542. 



Neutron sources 

Neutron sources were extensively used for detector calibra- 
tion, both to define an absolute energy scale and to mea- 
sure the neutron detection efficiency. A 252 Cf source was 
the most frequently used. This nucleus undergoes spon- 
taneous fission, simultaneously emitting prompt neutrons 
and 7-rays below 10 MeV. The neutron multiplicity obeys 
Poissonian statistics with an average value 3.787 [SHI and 
the kinetic energy is Maxwell-distributed, with an average 
energy equal to 1.3 MeV. 

We also used a Am/Be source. This emits neutrons, 
part of which follow the emission of a 4.4 MeV prompt 
7-ray (see jH] and references therein for details). 

With the source inserted in the centre of Region I, al- 
most all the neutrons were contained and captured: some 
(« 15%) were captured on Hydrogen, the majority on 
Gadolinium. The neutron binding energy was released in 
the form of 7-rays (one 2.2 MeV 7 for H, an average of 3 
7's and Y^ 7 E y w 8 MeV for Gd). In order to record events 
due to neutron capture on H as well, we defined a different 
trigger logic, by-passing the Llhi threshold condition; the 
L2 condition relied on the occurrence of two signals both 
satisfying the Lllo threshold. The "positron" role was, in 
this case, played by prompt 7's, while the second LI trig- 
ger was provided by the first neutron capture (if within 
the 100 /zs coincidence time- window) . 

The source could be inserted in the detector through 
two vertical calibration pipes: a central pipe, along the 
symmetry axis of the detector, and a parallel one, in Re- 
gion II, just at half distance between the geode and the 
target boundary, as shown in Fig. ED It was therefore 
possible to study the detector response anywhere along 
the z-axis in both regions. In particular, as we will see 
in the next section, calibration runs at different, known, 
positions allowed accurate tuning of the reconstruction al- 
gorithm and comparison with Monte Carlo forecasts. 



7 sources 

Daily calibration runs were performed by means of a 60 Co 
source, a well-known 7-emitter providing a low-energy cal- 
ibration point at 2.5 MeV (the sum of 1.17 and 1.33 MeV 
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Fig. 18. Location of laser flashers and calibration pipes in the 
detector. 



simultaneous 7-lines). As we will see in the following sec- 
tions, these runs were particularly useful to follow the de- 
tector evolution (a critical point because of the aging of 
the Gd- loaded scintillator), and to tune and check the trig- 
ger efficiency. 

The neutrino trigger was not suited for acquiring these 
events, since the delayed neutron signal was missing. In 

this case, the L2 trigger was made to coincide with the 5 Data acquisition 
Lllo trigger. With such a trigger we also recorded a great 
number of events due to 7-radioactivity natural background 



In this case the L2 was given by the internal laser trig- 
ger after a delayed of « 200 ns relative to the light emis- 
sion. 



4.4 The neural network trigger 

This trigger [S3] was intended to work as a fast acquisi- 
tion filter for selecting events inside a fiducial volume, and 
secondly, it could be used for triggering upon coincidence 
of two Lllo triggers, therefore acquiring those neutrino 
events followed by neutron capture on Hydrogen (instead 
of Gadolinium). 

This technique is not essential in the case of CHOOZ, 
since the amount of H-captures is w 15% of the total 
neutrino interactions in the target, but it could be im- 
portant in on-going reactor experiments (namely Kam- 
LAND), where only a high-flash liquid scintillator (not 
Gd-loaded) is used. 

The neural network algorithm was implemented on a 
VME parallel processor which reconstructed in w 180 (is 
the energy and position of each event from the informa- 
tion provided by the NNADC digitizing system (see next 
section). 

Due to the progressive deterioration of the Gd scintil- 
lator the potentialities of this trigger could not be fully 
exploited in this experiment. 



The data acquisition system consisted of various and dis- 
tinct systems controlled by a VME-based processor: 



Laser 

A laser based calibration system was also available. It con- 
sisted of a nitrogen UV Laser emitter (A = 337 nm) with 
a repetition rate ~ 10 Hz and a time jitter of less than 
1 ns between the trigger and the light emission. The light 
intensity could be selected by using two variable neutral 
density filter wheels rotated by remote-controlled stepping 
motors. The attenuated light was fed to any one of six lo- 
cations (remotely selectable) , using UV quartz fibers (only 
two out of the four Veto flasher positions are shown in 
Fig. 1181) . At each position the bare fiber was terminated in 
a quartz vial filled with a dilute p-terphenyl solution, com- 
prising the flasher. The fluorescence process shifted the 
laser wavelength to the same emission spectrum produced 
by scintillation, resulting in an isotropic point source hav- 
ing the same attenuation properties, time characteristics 
and PMT quantum efficiency as light generated by ener- 
getic particles. Time and pulse height calibrations of the 
detector could be carried out and monitored ona pulse- 
to-pulse basis. The monitors consisted of a photodiode, 
used to monitor laser intensity, and a separate fiber that 
directed light passing through the attenuator into a mon- 
itor PMT located next to the laser. Each of the monitors 
could be calibrated relative to an absolute standard. 



NNADC -» ADC units (VME), 

SWFD -> Slow (20 MHz) WFD units (VME), 

FWFD -> Fast (150 MHz) WFD units (FastBus), 

FBADC -» ADC units (FastBus), 

FBTDC -> TDC units (FastBus). 

The PMT signals were fanned out and fed into each digi- 
tizing system, according to the diagram shown in Fig. ^3 
Since the PMT's were used in the grounded cathode con- 
figuration, their signals were fed to the circuits via AC 
couplings. For the first three systems, the signals from 
groups of PMTs were fanned in; for the NNADC's and 
SWFD's the PMT signals were linearly added by groups 
of 8, a PMT patch. This reduced the number of channels 
to 24. For the FWFD's the PMT groups varied from 4 to 
8 during the experiment. 



5.1 The on-line system 

The acquisition system combined different bus standards 
(see Fig. 1201 : the digitizers listed above used both the 
VME and the FastBus standard, while the trigger elec- 
tronics (including discriminators and logic units) was CA- 
MAC based. The data were managed by a VME, OS-9 
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Fig. 19. Electronics layout for the CHOOZ experiment, in- 
cluding front-end, trigger and digitizing modules. 
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units 



__ 



discri Themi: 



Fig. 20. On-line system architecture, with special reference 
to the different bus standards (VME, FastBus, CAMAC) and 
their interconnection (VIC,VMV). 



operating processor, whose central unit consisted of a Mo- 
torola 68040 microprocessor mounted in a CES FIC 8234 
board. After completing the event readout, the data were 
sent through Ethernet to a dedicated SUN/Unix station 
and written on disk. A Lab VIEW run controller was the 
interface between the user and the SUN; the controller also 
provided real-time information about the run and data 
quality. 



- 1 input/output register Caen V513: 

- 1 Adaptive Solutions CNAPS parallel processor board. 

The scaler and flash- ADC provided redundant informa- 
tion on the relative timing of events. The flash- ADC board 
was also used to sample the trigger signals (Lllo, Llhi, L2 
and VETO). The other two boards were used to produce 
the Neural Network trigger. 

The SWFD's 



5.2 The digitizing electronics 
The NNADC's 

The signals from the 24 PMT patches were sent to two 
VME 12-bit ADC banks, each composed of three 8-channel 
Caen V465 boards. Two banks were necessary to avoid 
any dead time between positron and neutron pulses, the 
time needed for the charge integration being w 14 [is. The 
charge was integrated within a fast 200 ns gate. An ex- 
ternal logic alternately switched the gate to the banks on 
the occurrence of an Lllo trigger. The digitized ADC val- 
ues were arranged in a FIFO, thus allowing an internal 
multi-hit buffering in each bank. The ADC readout was 
performed by a Motorola 68030-based Themis TSVME133 
VME processor board (referred to as "Themis" in Fig. 120(1 . 
The TSVME133 was usually in a wait status until it re- 
ceived a VME interrupt generated when one of the two 
ADC banks was not empty. 

The NNADC circuitry also included: 

- 1 scaler Caen V560N: 

- 1 flash-ADC Caen V534: 



The Slow WaveForm Digitizer system used a set of 5 
flash- ADC Caen V534 boards, with 20 MHz internal clock 
speed. These digitizers recorded the 30 Geode plus Veto 
patch signals from 100 /is before to 100 /is after the L2 
trigger. All the input signals needed a proper stretching 
in order to fit to the inherent 50 ns resolution. Therefore 
these digitizers were essentially peak sensing devices and 
could not be used for pulse shape analysis. Extra channels 
recorded all the detector triggers. 

The FWFD's 

With their 155 MHz sample clock, these digitizers were 
the only equipment which could permit the study of pulse 
shapes. These units were equipped with an 8-page mem- 
ory, the depth of one page being 206 ns, in order to record 
up to 8 events before the stop. 

The FBADC's 

The FastBus ADC system had four 96-channel Lecroy 
1885F boards to separately measure the anode charge of 
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each geode PMT. Even in this case, the boards were ar- 
ranged in two alternating banks in order to record both 
events associated with the L2 trigger. 



The FBTDC's 

The FastBus TDC's recorded the time information of each 
geode PMT. Referring to Fig. ED the input signals con- 
sisted of the logical (NIM) pulses generated by the trigger 
discriminators at minimum threshold. This system con- 
sisted of two 96-channel Lecroy 1877 multi-hit units able 
to record up to 16 hits in one single channel with 2 ns 
resolution. 



6 Detector calibration 

The reconstruction methods were based, for each event, 
on the VME ADC's digitization of the signals from the 24 
patches which grouped the 192 PMTs looking at Region 
I. It was checked that this grouping does not significantly 
affect the quality of the reconstruction. 

A good determination of the position and energy of 
each event and of the detection efficiency depends to a 
large extent on the knowledge of the main parameters 
of the detector (scintillator light yield and attenuation 
length, PMT and electronics gains). 

6.1 PMT gain 

The stability of the PMT gain was checked by periodic 
measurements of the single photoelectron pulse height spec- 
trum of each PMT. Each spectrum was obtained by record- 
ing the PMT counting rate, out of a discriminator, as 
a function of its threshold. The data were taken by a 
Lab VIEW code running in parallel with the main ac- 
quisition. No artificial light source was needed, since the 
single counting rate, at the one-pe level, due to the 7- 
radioactivity and to the dark noise, was ~ 1 kHz. Ex- 
amples of these curves are shown in Fig. [21 The sin- 
gle pe spectra, obtained by differentiating the counting 
rate curves, were then fitted to get the PMT gains; val- 
ues around 30 mV for the geode PMTs and 5 mV for the 
Veto ones were needed to match our front-end electron- 
ics. Fig. EH shows the distribution of the average single 
pe pulse height for the whole of the geode PMTs and its 
evolution since the beginning of the data taking. One con- 
cludes that the PMTs stabilized about one month after 
burn-in and that the asymptotic gain was ~ 10% lower 
than the starting value. 



6.2 Determination of the photoelectron yield 

A similar method was used to determine the photoelectron 
yield for events at the detector centre. We used the laser 
flasher at the detector centre as a calibrated light source. 
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Fig. 21. Pulse height spectra for a sample of four PMTs. 
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Fig. 22. Distribution of the single pe peak for all PMTs (left) 
and its time evolution since the start of data taking (right). 



The laser provides its own trigger, so the counting rate is 
not affected by radioactivity. 

An estimate of the number of photoelectrons N l de- 
tected by each PMT comes from the probability of having 
zero photoelectrons when the average number is N l . From 
Poisson statistics, we obtained for the i-th PMT 

TV? 

iV l = -log(l-— (18) 

e l N shot 

where iVj^ t , iV s /> ot are respectively the counts over thresh- 
old for that PMT and the number of laser shots. (|T%)l 
includes also a detection efficiency s l , equal to the proba- 
bility for a single pe pulse height to be over threshold, 
which can be evaluated for each PMT from the pulse 
height spectra shown in Fig. [21 The light intensity was 
kept at a level of about 0.5 pe/PMT so as to increase the 
sensitivity of the method. This number was divided by 
the energy deposit corresponding to the light intensity in 
order to determine the absolute light yield. As an energy 
reference we used the 2.5 MeV 60 Co "sum" line. By av- 
eraging the number of photoelectrons all over the PMTs, 
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we obtained (0.65 ± 0.03) pe/MeV/PMT and a yield of 
(125 ±5) pe/MeV. 

6.3 ADC calibration in a single photoelectron regime 

Calibration runs using the laser at a low intensity were pe- 
riodically performed to test the single photoelectron gain 
and to calibrate the VME ADCs. The good ADC resolu- 
tion and the high PMT gain allow us to distinguish the 
contributions to the ADC spectra due to different numbers 
of photoelectrons. A sketch of these spectra is presented in 
Fig- 123 The fit function results from the sum of a pedestal 
and a number of pe-distributions; each of these terms is 
represented by a Gaussian function (weighted according 
to Poisson statistics) whose peak position is linear with 
the number of photoelectrons. 








ADC channels 



Fig. 23. ADC spectra obtained with a laser calibration, where 
the single pe peak is clearly visible. The fitting function results 
from the sum of the pedestal and the first three photoelectrons. 



Apart from the ADC calibration, the fit also provides 
an independent determination of the number of photo- 
electrons collected by each PMT patch which is consistent 
with the measured photoelectron yield. 



6.4 Light attenuation in the Gd-loaded scintillator 

The degradation of light transparency in Region I (see 
H3.3jl was regularly monitored in the apparatus by record- 
ing the number of photoelectrons associated with the 60 Co 
7-line, as shown in Fig. 1241 



.e 



Period I 



Period II 



Co peak (pedestal subtracted): t ~ 720 d 



* Threshold (pedestal subtracted) 



Days since Mar 12 1997 

Fig. 24. Peak associated with the 60 Co 2.5 MeV line, as a 
function of time, as measured by means of a Lecroy QVT. The 
detected charge follows an exponential decrease, with decay 
time « 720 d. 



The light absorption of Region I scintillator was peri- 
odically measured throughout the acquisition period. The 
method we used consisted in displacing a light source 
along the calibration pipe and in recording the charge 
detected by the top and bottom PMT patches. If we as- 
sume an exponential light attenuation (which is a good 
approximation for a light path longer than 10 cm, as al- 
ready remarked), these charge values are related by the 
simple expression 



Qb f?B AGd 



(19) 



£2t,b being the solid angle subtended by the top (bottom) 
patch PMTs and dr.s the average distance between these 
PMTs and the source position. 

For this measurement we used a radioactive source, 
namely 252 Cf. This source was preferred to a 7 source 
(such as 60 Co) since the double signature, provided by 
the prompt 7's and the neutron captures, makes the iden- 
tification of source events much easier. The 252 Cf neutron 
emission and capture, the involved light paths and solid 
angles were simulated by the MC method. 

The results obtained from different measurements are 
displayed in Fig. 1251 superimposed on the exponential 
best- fit curves according to (|19|l . The charge values in use 
are those corresponding to the 2.2 MeV 7-line due to the 
neutron capture on Hydrogen. We did not use the 8 MeV 
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line in order to reduce systematics due to ADC satura- 
tion effect which might have arisen when the source ap- 
proached the edges of Region I. The fitted Xcd values are 



active electronic components in the circuits treating the 
signals from the PMTs. Amplification balancing factors 
were obtained three times a week using the ADC charge 
values corresponding to the 8 MeV line peak generated 
with the 252 Cf source at the detector centre. 
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Fig. 25. Attenuation length measurements at different stages 
of the data taking period. 



plotted versus time in Fig. EEl The time evolution of the 
absorption length was fitted by the empirical function 11611 
(already used in the laboratory test). By taking the Xcd 
values, at the beginning and at the end of the experiment, 
one can estimate a 35% reduction in the photoelectron 
yield, in agreement with the measurement of the QSUM 
for events at the detector centre (Fig. 124ft . 



7 Event reconstruction techniques 

7.1 The standard minimization algorithm 

The standard algorithm uses a maximum likelihood method 
to reconstruct the energy E and the vertex ~>t of an event. 
The likelihood is defined as the joint Poissonian probabil- 
ity of observing a measured distribution of photoclectrons 
over the 24 patches for given (E, ~x) coordinates in the 
detector. So, for an event occurring at time t after the 
start of data taking, we can build a likelihood function as 
follows: 
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C(N; N) = H P(Nj;Nj(E, 1?, t)) = J] 



N Nl - 
N 3 \ G 



(20) 

where Nj is the observed number of photoelectrons and 
Nj the expected one for the j-th patch, given an event 
(E,~x,t). The reason for using a Poissonian instead of 
Gaussian statistics is due to the frequent occurrence of low 
energy events, with low number of photoclectrons detected 
by some PMT patches. 

The values Nj are obtained from the data recorded by 
the VME ADC's by applying the reference ADC gain go 
and the balancing factors fj : so 




Ut) = y(l + at) 

\ = (587 ± 33) cm 
a = (4.2 ± 0.4) 10 3 d 1 
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Fig. 26. Xcd versus time with best-fit function superimposed. 



6.5 Electronics amplification balancing 

The electronic gain may differ from patch to patch and 
slightly vary with time because of the behaviour of the 



9ofj(t) 



(21) 



The predicted number of photoelectrons for patch j is 
computed by considering a local deposit of energy, result- 
ing in a number of visible photons which are tracked to 
each PMT through the different attenuating Region 1 and 
2 scintillators. Therefore 



Nj = aEr)J2 
where 



f2 jk (~>t) 



fe=i 



-in 



exp 



rfijfc(~x) _ d 2 jk(~x) 
Xad(t) Xffi 



(22) 



E is the energy deposited in the scintillators, 

a is the light yield of the scintillator, 

rj is the average PMT quantum efficiency, 

Qjk is the solid angle subtended by the k-th PMT 
from the position, 

dijk is the path length in Region I, 

dijk is the path length in Region II, 

Xcd is the attenuation length in Region I scintillator, 

\ui is the attenuation length in Region II scintillator. 
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To reduce computing time PMTs are considered to be 
flat and the solid angle is approximated by the following 
expression 



n jk = 2tt 1 1 



V k 



d %+ r lMT C0S < 



(23) 



rpMT being the PMT photocathode radius, 9 the angle be- 
tween the event-PMT direction and the inward unit vector 
normal to the PMT surface and djk = rfyfc + dijk- 

Instead of directly using Eq. i|20[) it is more convenient 
to exploit the "likelihood ratio test" theorem to convert 
the likelihood function into a form which obeys the x 2 
distribution[SHI. We let Nj be the best estimate of the 
true (unknown) photoelectron distribution and form the 
likelihood ratio A defined by 



A = 



£(N;N) 
C(N;N) 



(24) 



The "likelihood ratio test" theorem states that the "Pois- 
sonian" x 2 ; defined by 
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X 2 = -21ogA = 2]T[TV 

3=1 



N- 

, -JV 3 -+^log(^)], (25) 



asymptotically obeys a chi-square distribution jHZj. It is 
easy to prove that the minimization of \ 2 is equivalent 
to maximization of the likelihood function, so that the 
X 2 statistic may be useful both for estimating the event 
characteristics and for goodness-of-fit testing. 

We used the MIGRAD minimizer provided by the 
CERN MINUIT package [£E] to minimize The search 
for the minimum x 2 proceeds through the computation of 
the first derivatives of l|25(l . This routine proved very pow- 
erful, provided the starting values for the fit parameter 
are accurate. We studied the x 2 profile by reconstruct- 
ing Monte Carlo generated events. Several relative min- 
ima were found, most of them differing by more than la 
from the generated (E, ~x") coordinates. This is the reason 
why a suitable choice of these starting values is crucial 
in event reconstruction. In our case, an indication of the 
event position comes from the asymmetry of the charge 
distribution: no significant asymmetry will be visible for 
events at the centre while the charge distribution will be 
more and more asymmetric for events approaching the de- 
tector boundary. We subdivided the PMT patches into 6 
"superpatches" , two patches for each coordinate, and built 
the starting point for the i-th coordinate according to the 
following equation: 




1,2,3, 



(26) 



where the indices H — refer to the opposite superpatches 
of the i-th axis and D l is the half size of the detector 
along that axis. Once the corresponding to the start- 
ing position is known, the starting energy value is obtained 



from ill' 21) after replacing ~x with 'x'o and Nj with Nj. 
The event reconstruction was tested by analysing events 
generated with calibration sources in various positions in- 
side the detector. Let us review the results. Distributions 
of reconstructed events generated by the laser flasher at 
the detector centre are presented in Fig. [23 The stan- 
dard deviation of the distributions shown gives an in- 
dication of the resolution, both in energy and position. 
The fit yields a x ss 4 cm for each coordinate and an en- 
ergy resolution oe — 0.33 MeV at an equivalent energy of 
8.1 MeV in which the statistical fluctuations in the num- 
ber of photoelectrons prevail. Photoelectron statistics also 
affect the position resolution. The effect is clearly visible 
in the data taken with the same laser flasher for three 
different light intensities (in Fig. [23 a larger filter atten- 
uation corresponds to a lower intensity). As can be seen 




Fig. 27. Comparison of position and energy distributions for 
runs with the laser flasher at the detector centre, corresponding 
to three different light intensities. 



in Fig. |23 the average z coordinate is displaced by 7 cm 
from the nominal position; this is due to the effect of the 
shadow of the flasher already mentioned; in this case, the 
flasher points downward, thus producing a shadow in the 
upward direction and displacing the event below the true 
source position. In Figs. |2H1 12H1 we report the results of 
the calibration runs with the 252 Cf source at two differ- 
ent positions [z = 0, —80 cm) along the calibration pipe. 
The position distribution is Gaussian for all the coordi- 
nates, with a x w 19 cm. An equivalent number of neutrons 
was generated, at each position, by our Monte Carlo code 
and similarly analysed. The figures also display the Monte 
Carlo data to emphasize the agreement between data and 
expectations. 

The energy spectra show the energy lines due to neu- 
tron capture on Gd and H. The shape of the Gd-capture 
energy spectrum, as pointed out in Fig. 1301 results from 
the superimposition of 7 lines due to neutron capture on 
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Fig. 28. Data and Monte Carlo distributions of neutron events 
with the 252 Cf source at the detector centre. Dashed and solid 
lines in the bottom figures represent respectively the top and 
the side view of the edge of the target and the geode. 



157 Gd and 155 Gd (see Tab.EJfor a reference). Due to the 
scintillator saturation (see ^I3.5[) the fitted peak values are 
~ 2.5% lower than the nominal ones. 

7.2 Reconstruction problems 

We have just seen that event reconstruction gives good 
results, for calibration runs, as long as the source position 
is inside Region I. Unfortunately this is not always true 
and more and more problems arise for events closer and 
closer to the geode surface. 

The main problem concerns the 1/r 2 divergence of the 
light collected by one PMT; the exponential light atten- 
uation, entering formulae I|19I22[1 . becomes inaccurate in 
the vicinity of the PMTs. In such a case, the approxima- 
tion of a flat PMT surface is also no longer adequate to 
evaluate the solid angle by means of (|23l) . These effects 
are shown in Fig. 1311 The neutron capture energy (both 
for Gadolinium and Hydrogen) is overestimated at recon- 
structed distances smaller than ~ 30 cm from the geode 
surface. 



Fig. 29. Same as before, with the source at z — —80 cm. 



Fig. also shows one further weakness inherent in 
the minimization procedure. As remarked above, the MI- 
GRAD minimizer heavily depends on the knowledge of the 
first derivatives of (|25|l with respect to the fit parameters 
and fails if this computation is not accurate enough. This 
situation arises for events near the acrylic vessel, where 
the first derivatives of l|25|l exhibit a discontinuity due to 
the different light attenuation of the Region I and Region 
II scintillators; as a result, a discontinuity in the energy vs. 
distance profile is found around the vessel surface position. 
A slight improvement was obtained by calling the CERN 
SIMPLEX routine when MIGRAD failed. This minimizer, 
although much slower and in general less reliable than MI- 
GRAD, is more suitable in this case since it does not use 
first derivatives and it is not so sensitive to large fluctua- 
tions in the function value. 

Fig. 1321 shows the distribution of 252 Cf neutron events 
generated at z = —120 cm. The discontinuity in the z dis- 
tribution is also present for Monte Carlo generated events, 
thus proving that this is an inherent feature of the mini- 
mization procedure, and was properly taken into account 
in evaluating the neutrino detection efficiency. 
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Cf source, z = -120 cm 
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Fig. 30. Reconstructed energy spectrum for events associated 
with neutron capture on Gd. Contributions from 7-lines at 
7.94 MeV (capture on 157 Gd) and 8.54 MeV (capture on 155 Gd) 
are singled out. The double-Gaussian fit parameters are also 
shown. 



„ 20 
> 



" 18 

S? 
1 

u 16 

e; 

o 

a 

g 14 



aaylic 
vessel '■ 




20 40 60 80 100 120 140 160 
distance from edge (cm) 



Fig. 31. Energy versus distance from edge for neutron events 
generated by the 252 Cf source along the calibration pipe. The 
energy is flat whenever the reconstructed position is more than 
30 cm from PMTs. 



8 Neutrino event selection and background 
rejection 

Neutrino events must be identified in a diverse and much 
more numerous background. Off-line selections, based on 
the reconstructed event variables, were chosen to enhance 
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Fig. 32. Distributions of neutron events with the 252 Cf source 
at z = -120 cm. The discontinuity in the z distribution at the 
vessel surface is visible also in Monte Carlo generated events. 



the signal/background ratio; at full reactor power this ra- 
tio is > 20. We performed a careful evaluation of the effi- 
ciencies associated with the selection criteria. The analysis 
procedure is reviewed and the characteristics of the neu- 
trino event sample are presented. 



8.1 The data sample 

The data acquisition period extended over 450 days, from 
12 March 1997 till 20 July 1998; about 2000 runs were 
taken, including standard neutrino runs (about 1/4 of the 
total) and daily calibrations. The experiment stopped tak- 
ing data five months after the last reactor shut-down (8 
February 1998) when it became clear that, due to prob- 
lems related with the cooling system, neither reactor would 
resume its normal operating conditions for at least one 
year. 

Data taking is summarized in Tab. ^ the total ther- 
mal energy released (which can be intended as a neutrino 
integrated luminosity) is also listed. The total live time 
amounts to sa 340 d, 40% of which was with both reac- 
tors off. The power evolution of the CHOOZ reactors is 
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illustrated in Fig. 1231 The set of power values for both 
reactors almost continuously covers the entire range up to 
full power. It is worth noting that this is a very unique 
feature of the CHOOZ experiment. 



i 




Fig. 33. Power (top) and burn-up (bottom) evolution for 
CHOOZ reactors. Both have been off since February 1998. 



Only one of the two CHOOZ reactors was on for at 
least 80% of the total live time of the experiment. This 
allowed us to extract the separate contribution of each 
reactor to the neutrino signal and to perform, thanks to 
their different distances from the detector, a two-distance 
oscillation test. 



8.2 Candidate event selection 

Radioactivity background events are greatly reduced by 
applying loose energy cuts to the neutron-like signals; 
events are selected if QSUM n > 13000 ADC counts, which 
roughly corresponds to a 4 MeV energy deposit at the de- 
tector centre. The residual events (« 7.2 • 10 5 over a total 
number of « 1.2 • 10 7 L2 triggers) are then reconstructed 
by the standard minimization procedure described above. 

An analysis of this preliminary sample gives a first il- 
lustration of the properties of the neutrino signal (reactor- 
on data) and the associated background (reactor-off data). 
Figures |23 show the correlation of neutron-like vs. 
positron-like energy for this sample. Neutrino events, as 
indicated in Fig. EH populate a region in the (E B +,E n ) 
plot delimited by E c + < 8 MeV and 6 < E D < 12 MeV. 
Background events, depending on their position in these 
scatter plots, are classified in the following categories: 

A) events with E c + < 8 MeV and E a > 12 MeV: 
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Fig. 34. Neutron versus positron energy for neutrino-like 
events collected during the reactor-on period. A preliminary 
cut to the neutron QSUM is applied to reject most of the ra- 
dioactivity background. 
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Fig. 35. Same as before, with neutrino-like events collected 
during the shut-down of both reactors. 



B) events with E c + > 8 MeV and E D > 12 MeV: 

C) events with E c+ > 8 MeV and 6 < E n < 12 MeV: 

D) events with E c+ < 8 MeV and E n < 6 MeV: 

The neutron energy distribution of C) events presents 
a clear 8 MeV peak, typical of the neutron capture on 
Gadolinium, still persisting in the reactor-off data. These 
events can then be interpreted as a correlated background 
associated with high energy spallation neutrons from cos- 
mic ray interactions in the rock surrounding the detec- 
tor; the neutrons entering the detector are slowed down 
to thermal velocities via elastic scattering on protons and 
then captured; the proton recoil signal, whose energy spec- 
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trum is basically flat and extends to high energies, mimics 
the positron signal. This interpretation is confirmed by 
the neutron delay distribution. An example is shown in 
Fig. 1361 where events in sample C) follow an exponential 
decay distribution whose life time r = (30.5 ± 1.0) /its is 
correct for neutron capture in the Gd-doped scintillator. 



overburden. However, the neutron signals associated with 
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Fig. 36. Distribution of positron-neutron delay for the differ- 
ent event categories. The best fit curves are also drawn and 
the relative parameters indicated. 
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Fig. 37. Neutron versus positron energy for neutrino-like 
events selected from the preliminary sample by applying the 
"topological" cuts here indicated. 



this background component are an important tool to fol- 
low the energy calibration stability throughout the exper- 
iment. The results are shown in Fig. 1381 where the av- 
erage neutron energy, obtained by a Gaussian fit of the 
Gadolinium capture peak, is plotted vs. the run number. 
The relative stability throughout the data taking period 
(drift at a 0.8% level) is an independent verification of the 
reliability of the adopted reconstruction technique. 



Events in categories A) and D) show a flat delay dis- 
tribution, therefore implying an accidental coincidence of 
uncorrelated signals. In both cases, the positron signal is 
faked by a low energy (E 1 < 3MeV) radioactivity event. 
The neutron signal may be associated with either another 
radioactivity event, as in category D), or a high activity 
signal (most likely due to a proton recoil) as in category 
A). 

Finally, the delay of B) events is exponentially dis- 
tributed with a lifetime r = (2.8 ± 0.4) /is which is com- 
patible with that of a muon decay at rest. These events 
can be associated with residual cosmic muons stopping in 
the detector and then decaying. Both the muon energy 
loss and the Michel electron energy are much higher than 
the typical energy of a reactor neutrino interaction; these 
events can be rejected by applying an energy selection 
both to the positron-like and the neutron-like signals. 

The accidental background can be significantly reduced 
by applying fiducial volume cuts. Less than 10% of the 
events in regions A) and D), as shown in Fig. 1371 survive 
the selection cuts (distances from the geode d > 30 cm 
for both positron and neutron and relative distance < 
100 cm). Conversely the correlated background is more 
difficult to eliminate since these events exhibit the same 
features as neutrino interactions. This is the reason why 
the final background rate is dominated by the correlated 
component, in spite of the shielding provided by the rock 



<E> = 7.64 MeV 
ov = 0.50 MeV 




Fig. 38. Neutron energy distribution for correlated back- 
ground events (left) and average E n vs. run number (right). 



8.3 Final selection 

Both the energy and the topological cuts were studied 
and optimized by relying on the Monte Carlo simulation 
of neutrino events. These predictions were cross-checked 
with the calibration data to gain confidence in the final 
efficiencies. We adopted the following criteria for selecting 
neutrino events: 

1) Positron energy: E c + < 8 MeV: 

2) Neutron energy: 6 < E n < 12 MeV: 
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3) Distance from geode boundary: d c + > 30 cm, d n > 
30 cm: 

4) Relative positron-neutron distance: d c + n < 100 cm: 

5) Neutron delay: 2 < At c+n < 100 ^s: 

6) Neutron multiplicity: N n = 1. 

The cut on the positron-like energy is chosen to accept 
all possible positron triggers. The Lllo trigger threshold 
limits the lower end while the upper limit was set to 8 MeV 
since the probability of having a larger positron energy is 
negligible (< 0.05%). 

Cut 2) selects the neutrino events associated with the 
neutron captures on Gadolinium, which are more than 
80% of the total. The lower cut introduces an additional 
inefficiency due to 7 rays escaping the containment re- 
gion. As illustrated in Fig. 6 MeV is a suitable choice 
to separate neutrino events from the residual low energy 
uncorrelated background, thus optimizing the signal to 
noise ratio. 

The cut on the neutron delay covers about three cap- 
ture times for neutrons in Region I. The 2 fjs lower cut was 
introduced to reduce the effects of the signal overshoot in- 
herent the AC coupling of PMT bases and front-end elec- 
tronics. These effects are particularly troublesome in the 
case of NNADC's, which integrate current signals of both 
polarities. 

The correlated background can be reduced by applying 
a cut on the secondary particle multiplicity 2 . Muon spalla- 
tion processes usually generate several neutrons, therefore 
more than one particle is likely to enter the detector and 
give a detectable signal. 

Events satisfying the selection criteria will be referred 
to as neutrino event candidates from now on. The next 
section reviews the efficiency for selecting neutrino inter- 
actions and the background rejection resulting from indi- 
vidual cuts. 



8.4 Positron efficiency 
Energy cut 

The expected positron energy spectrum is obtained by 
folding the kinetic energy spectrum coming from @ with 
the detector response function. Two problems then arise 
in evaluating the positron efficiency. Firstly, the neutrino 
spectral shape slightly varies along a reactor cycle as a 
consequence of the fuel burn-up; however, the informa- 
tion daily provided by E.D.F. allows us to accurately fol- 
low this variation. Secondly, the energy cut operated by 
the Lllo trigger increases with time because of the aging 
of Region I scintillator. A daily check of the equivalent en- 
ergy threshold is required to account for this effect. The 
procedure we followed is shown in Fig. A Lecroy qVt 
recorded the energy spectra due to a 60 Co source placed 
at the detector centre. Two spectra were taken, the first 
one with an external trigger provided by the Lllo signal 

2 The positron (or primary) signal is associated with the LI 
trigger preceding the L2. Every signal associated with or oc- 
curring after the L2 is referred to as a secondary event. 



and the second with an internal trigger set at a lower 
threshold. In addition to the 2.5 MeV "sum" line, each 
spectrum exhibits a low energy tail due to the energy loss 
in the calibration pipe. The ratio of the two spectra (after 
background subtraction) yields the Lllo trigger efficiency 
curve. The equivalent energy threshold results from fitting 
this curve by an integral Gaussian function. 
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Fig. 39. Determination of the Lllo equivalent energy thresh- 
old. The top figures show the QSUM spectra measured with 
a 60 Co at the detector centre by means of an internal and ex- 
ternal triggered qVt (the corresponding background is super- 
imposed). The central plot shows the background subtracted 
spectra. The bottom histogram, displaying the efficiency curve, 
follows an integral Gaussian function whose parameters are in- 
dicated. 



The energy threshold varies with time as illustrated in 
Fig. UJ]| Its behaviour can be described by a linear func- 
tion. Discontinuities arise in coincidence with the thresh- 
old resettings, which were needed to bring the Lllo thresh- 
old back to the optimal value. 

More extensive 60 Co calibrations were periodically per- 
formed at different positions of the source along the cen- 
tral pipe, in order to measure the z dependence of the 
energy threshold. The results are shown in Fig. |2] at 
four different stages of data acquisition. The overall de- 
pendence on the position of the energy threshold can also 
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Fig. 40. Equivalent energy threshold at detector centre as a 
function of time. The jumps visible here are due to the thresh- 
old setting. 



be predicted by a Monte Carlo simulation of the QSUM 
and NSUM behaviour. The attenuation length values are 
the same as those used in the event reconstruction. The 
threshold measurement always fits well with expectations, 
as evident from the same figure. We can therefore rely on 
Monte Carlo predictions to determine the energy thresh- 
old at each position in the detector, where the source can- 
not be placed. An expression for the positron threshold as 
a function of time and position can be derived by com- 
bining the daily threshold measurements at the detector 
centre with Monte Carlo evaluations. An explicit form of 
this threshold is given in cylindrical coordinates (/?, z, 4>) 
by the following equation: 

E thr (t, = E thr {t, 7 = 0) x (1 - a{t) P 2 ) x (1 - (3{t)z 2 ) 

(27) 

where the energy threshold at centre E t hr(t, x = 0) is 
extracted by interpolating the daily measurements shown 
in Fig. flUI and a(t),/3(t) are linear functions obtained by 
Monte Carlo. 



Distance from the geode cut 



The positron loss due to the selection 3) (distance from 
the geode boundary) was evaluated by a Monte Carlo sim- 
ulation of neutrino interactions. 10000 events were uni- 
formly generated in a volume including a 10 cm wide shell 
surrounding the target, so as to take into account spill- 
in/spill-out effects. The estimated efficiency proves to be 
almost independent of the scintillator degradation; we can 
then consider it as a constant and get the average value 
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Fig. 41. Energy threshold as a function of z for different val- 
ues of the attenuation length for the Gd-doped scintillator. 
The measurements obtained with the 60 Co source follow the 
expected behaviour. 



8.5 Neutron efficiency 
Neutron capture cut 

The neutron capture efficiency is defined as the ratio of 
neutrons captured by Gadolinium nuclei to the total num- 
ber of captures. This ratio enters the global neutron effi- 
ciency since the neutron energy selection excludes events 
associated with neutron captures on Hydrogen. This cap- 
ture efficiency was studied by means of a normal and a 
special tagged 252 Cf source, providing a fission tagging 
signal. The calibration pipe was removed for these partic- 
ular studies in order to avoid effects due to the iron con- 
tent in the pipe. The source was at the end of a plumb- 
line long enough to reach the bottom of the target ves- 
sel. Fission neutrons have a longer mean pathlength than 
neutrons from neutrino interactions (w 20 cm in scintilla- 
tor instead of ~ 6 cm) , so calibration data must be cou- 
pled with Monte Carlo predictions to correctly take into 
account the effects associated with the neutron spill-out 
(which are relevant for events approaching the boundary 
of Region I) . 
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Neutron events are associated with Gd captures if 4 < 
E n < 12 MeV. The ratio obtained includes a 1% correction 
due to Gd-capture events with visible energy below 4 MeV. 
The combination of data and MC studies results in an 
average efficiency over the entire detector of 



e Gd = (84.6 ±0.85)% 



(29) 



Energy containment cut 



The 6 MeV lower energy cut introduces an inefficiency in 
neutron detection due to the 7 rays from Gd-capture (by 
definition inside Region I) escaping the visible volume. 
Since this efficiency is almost independent of the neutron 
pathlength we can rely on the calibration data to evaluate 
this number; Monte Carlo predictions are used to cross- 
check the data. Since this effect is expected to be more 
relevant for events near the edge of Region I, a fine source 
scanning was made in this region. A special set of runs, 
was made with 2cm steps, between f.5cm and If. 5 cm, 
from the bottom edge of the acrylic vessel. 

The escape fraction was defined as the ratio of the 
number of events with 4 < E n < 6 MeV over those with 
4 < E n < f2MeV. The values turn out to range from 
2.9% at the centre to up to 6.6% at the bottom edge. By 
averaging over the target volume we get 
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Fig. 42. Neutron delay distribution measured with the Am/Be 
source at the detector centre (left) and at the bottom edge of 
the acrylic vessel (right). The time origin is denned by the 
4.4 MeV 7-ray. No time analysis cut is applied yet. 



Given the chemical composition of the scintillator, it is 
possible to rely on the Monte Carlo simulation to predict 
the efficiency related to the neutron delay cut in the case 
of neutrino induced events. The estimated loss due to the 
2 /its analysis cut amounts to f.6 ± 0.2%; the fraction of 
neutrons with At > f 00 /is is 4.7 ±0.3%. We therefore end 
up with: 

e At = (93.7 ±0.4)% (32) 



e En = (94.6 ±0.4)% 



(30) 



Distance from the geode cut 

We followed the same procedure used to evaluate the cor- 
responding positron efficiency. We obtain 

e da = (99.5 ±0.1)% (31) 

This value is not significantly affected by the scintillator 
degradation. 



8.6 Relative distance cut efficiency 

Again we used the reconstruction of 10000 Monte Carlo 
generated events to predict the efficiency due to the positron- 
neutron distance cut. We obtained 

£d ^ +n = (98.4 ±0.3)% (33) 

Again this efficiency is nearly independent of the time evo- 
lution of the Gd-doped scintillator. 



Delay cut 

The neutron delay cut was studied by using both the 252 Cf 
and an Am/Be source; the latter emits single neutrons 
only, preventing biases due to the ADC integration dead 
time arising at high neutron multiplicities. Fig. 1421 shows 
the delay distributions obtained with the Am/Be source 
placed at the centre and at the bottom edge of the tar- 
get. The neutron capture time at the centre can be used 
to determine the Gd content in the Region 1 scintillator 
(to be used in the Monte Carlo code). The fitted decay 
time t = (30.7 ± 0.5) fj,s corresponds to a mass fraction of 
(0.0940 ± 0.0015)%. With reference to Fig. Upright), at 
the bottom of the detector the neutron capture on Hydro- 
gen becomes more important and the fitted decay time 
increases accordingly. Both spectra show a low first bin 
content. This is due to neutron moderation time which 
corresponds to the relatively fast neutrons emitted by the 
Am/Be source. 



8.7 Neutron multiplicity 

The neutron multiplicity cut rejects neutrino events if a 
"spurious" Lllo trigger occurs along with the positron- 
neutron pair. Lllo triggers are mainly associated with 7- 
background events, « 97% of which have energies lower 
than the high threshold; so, the background above the 
Llhi threshold can be neglected in what follows. The error 
associated with this approximation is negligible. 
Let us consider all possible sequences: 

1) e + — n — 7 with t~ — t a < 100 fis: 

2) e+ - 7 - n: 

3) 7 — c + — n with t c + — t 7 < 100 fis: 
Case 1 

The L2 triggers on the (e + ,n) pair and the 7 signal is 
mistaken for a second neutron; the neutrino event is then 
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rejected because of cut 6). Given the 7 rate i? 7 , the prob- 
ability of such a sequence is 



1 — £1 — R~ t /\t v 



(34) 



(where At v is the 98 fis acceptance window) from which 
£1 is extracted. The average efficiency is (98.6 ± 0.3)%. 



Case 2 



If the positron energy exceeds the high threshold, the L2 
triggers on the (e + ,7) pair; then two neutron-like signals 
are detected (the "true" neutron plus the 7) and the neu- 
trino event is rejected by cut 6). If the positron energy is 
lower then the Llhi threshold, the L2 triggers on the (7, n) 
and the 7 is mistaken for a positron-like event. Since the 
7 signal is not correlated with the neutron, the probabil- 
ity of surviving the topological cuts 3) and 4) is quite low 
(< 2%) and can be neglected. As a conclusion, neutrino 
events occurring in such a sequence are always rejected, 
whatever the positron energy is. 

Let P n (t c ) be the neutron delay distribution shown in 
Fig. 021 an d t the time between the positron and the 7 
signals. The neutrino efficiency can be written as follows: 



£"2 



Pn(t c )dt c 



Rj dt 



R~yt c 



(35) 



where t c = (30.5 ± 0.5) /is is the average neutron delay. 



Case 3 



Like the above case, if the positron signal met the high 
trigger condition, L2 would trigger on the (7, e + ) pair and 
the selection criterium 6) would reject the event; otherwise 
the L2 would trigger on the (e + , n) pair and the neutrino 
events would be accepted. So 



Table 6. Summary of the neutrino detection efficiencies. 
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Fig. 43. Neutron energy spectra for reactor-on and reactor-off 
periods (left) and background subtracted spectrum compared 
to Monte Carlo expectations (right). 



£3 — 1 — RyfiiAt,, 



(36) 



fu = (0.45 ± 0.05) being the positron fraction above the 
high energy threshold. 



Combined efficiency 



The neutron multiplicity cut efficiency is obtained by mul- 
tiplying the three above efficiencies. We evaluated it on a 
run by run basis since i? 7 depends on the Lllo threshold. 
The average value is ei n = (97.4 ± 0.5)%. 

A summary of all the efficiencies is presented in Tab.|H] 
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Fig. 44. Same as before, for the positron-neutron distance. 

9 The neutrino signal 



Figures l4*3l through l45l present the final neutrino candidate 
sample, with all the selection cuts applied. A total number 
of 2991 neutrino candidates was collected, 287 of which oc- 
curred during the reactor-off periods. To properly compare 
these with expectations, the entire acquisition cycle was 
divided according to the dates of the threshold resetting, 
when variations in the background rate were expected. For 
each resulting period, the reactor-off background was nor- 
malized to the same livetime as the reactor-on spectra. 



9.1 Measured positron spectrum 

Fig. 1461 shows the complete measured spectra (reactor-on, 
reactor-off), obtained by summing the spectra collected 
during runs relative to different off-line periods; the re- 
sulting positron spectrum (reactor-on minus reactor-off) 
is presented in Fig. El The chosen bin width (0.4 MeV) is 
roughly adapted both to the statistics and to the energy 
resolution. 
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Fig. 45. Neutron delay distributions for reactor-on and 
reactor-off periods (left) and background subtracted spectrum 
compared to MC predictions (right). 
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Fig. 46. Experimental positron spectra for reactor-on and 
reactor-periods after application of all selection criteria. The 
errors shown are statistical. 



9.2 Predicted positron spectrum 

The expected visible energy positron spectrum at the de- 
tector position, for a mean reactor-detector distance L k , 
is given by: 

S k (E,L k ,e,5m 2 ) = -r^Np f h{L,L k )dL 

k J 

j <t(E c+ )S„{E„)P{E„, L, 6, 8m 2 )r{E c+ , E)e{E c+ ) AE C+ , 

(37) 

where 



E v , E e + 
N p 

a(E c+ ) 
S V {E V ) 
h{L,L k ) 

r{E c+ ,E) 

P(E„,L, 
e(E c+ ) 



?, 8m 2 ) 



are related by ® , 

is the total number of target protons 

in the Region I scintillator, 

is the detection cross section (fTTTf) . 

is the v c spectrum, 

is the spatial distribution function for 

the finite core and detector sizes, 

is the detector response function 

providing the visible e + energy, 

is the two-flavour survival probability, 

is the combined detection efficiency. 



We first computed the positron spectrum in the ab- 
sence of neutrino oscillations, by using the Monte Carlo 
code to simulate both reactors and the detector. The com- 
posite antineutrino spectrum was generated for each of the 
205 fuel elements of each reactor core; each antineutrino 
was assigned a weight according to the prescriptions given 
in H2.9I The interaction points were randomly chosen in 
the target; the positron and resulting annihilation photons 
were tracked in the detector, and scintillator saturation 
effects were included to correctly evaluate the positron 
visible energy. The resulting spectrum, summed over the 
two reactors, is superimposed on the one measured in 
Fig-E|to emphasize the fact that the data fit with the no- 
oscillation hypothesis; the Kolmogorov-Smirnov test for 
compatibility of the two distributions gives a 82% proba- 
bility. The measured vs. expected ratio, averaged over the 
energy spectrum (also presented in Fig. I47|l is 



R = 1 .01 ± 2.8%(stat) ± 2.7%(syst) 



(38) 



In the next section we will discuss how this ratio constrains 
the oscillation hypothesis. 



9.3 Background 

One of the main features of the experiment, the low back- 
ground, represents a difficulty in that it is hard to measure. 
Apart from the low statistics, another difficulty is that the 
background rate depends on the trigger conditions which, 
as discussed above, changed with time due to scintilla- 
tor aging and the positron threshold adjustments. Sepa- 
rate estimates of the background are then needed for each 
data taking period. Only 34 events were collected during 
the 1997 runs (periods 1 & 2, until 12 January 1998, date 
of the last threshold resetting: see Tab. below) with re- 
actors off, the total live time being 577 h; this implies a 
background rate of 1.41 ± 0.24 events per day. Most of the 
reactor-off statistics, amounting to 253 events, was col- 
lected during the 1998 run when the cumulative reactor 
power was very low; the reactor-off live time is 2737 h. 
The resulting background rate is 2.22 ± 0.14, about twice 
as large as in the 1997 run. 

The explanation for this variation relies on the low- 
ering of the NSUM threshold associated with the Lllo 
trigger. A lower number of hit PMTs implies a larger fidu- 
cial volume extended towards the PMT boundary, where 
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9.4 Correlated background 
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Fig. 47. (above) Expected positron spectrum for the case of 
no oscillations, superimposed on the measured positron ob- 
tained from the subtraction of reactor-on and reactor-off spec- 
tra; (below) measured over expected ratio. The errors shown 
are statistical. 



the event rate is dominated by the natural radioactiv- 
ity. Moreover, we saw that the reconstruction algorithm 
overestimates the energy of events near the PMTs (see 
Fig. Enj. Therefore a larger fraction of radioactivity events 
is shifted into the neutrino candidate event window shown 
in Figs. 021123123 As a result, the accidental component of 
the background is greatly enhanced. On the other hand, 
the correlated background is virtually unaffected by the 
low trigger conditions, since the neutron signal is much 
higher than the threshold. We verified that the correlated 
background did not significantly change throughout the 
experiment. 

As a further cross-check of the reactor-off estimates, 
we measured the background by extrapolating to zero the 
candidate yield versus reactor power. This will be the sub- 
ject of 103 



Fig. 1481 shows the energy distribution of e + -like signals as- 
sociated with the correlated background; these events are 
selected by applying all the criteria for neutrino candidate 
selection, except e + -like energy. Apart from a low-energy 
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Fig. 48. Energy distribution of e + -like signals associated with 
the correlated background. 



peak (due to the neutrino events and to the accidental 
background) , the spectrum has a roughly flat distribution 
extending beyond 30 MeV, with a slight increase at higher 
energies caused by the NNADC saturation. A Monte Carlo 
simulation of this spectrum is very difficult, since no reli- 
able transport code is available for spallation neutrons in 
the 10 -J- 100 MeV kinetic energy range. Moreover, the ob- 
served spectrum is affected by the scintillator saturation 
(which is relevant at low recoil proton energies). 

To verify the stability of the correlated background 
rate, we divided the complete run into three different pe- 
riods (see Tab.[7]below) corresponding to the dates of the 
threshold resetting. The spectra for each period were fit- 
ted in the 10 < E < 30 MeV range. We divided them by 
the live times (also listed in Tab. 0) of the three periods, 
and finally found: 



B, 



(0.156 ± 0.01) MeV" 1 d^ 1 for the 1st period, 
(0.158 ± 0.01) Me^ 1 d _1 for the 2nd period, 
(0.151 ± 0.01) Me^ 1 d _1 for the 3rd period, 

(39) 

confirming the stability of the correlated background rate 
throughout the experiment. 

A rough evaluation of the correlated background can 
be obtained by extrapolating the rate in the positron win- 
dow (E t hr < E < 8 MeV). Taking the average value from 
(I39|l and assuming Ethr = 1.5 MeV, we obtained a rate of 
(1.01 ± OM(stat) ± 0.1(syst))eventsd _1 . 
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9.5 Accidental background 

The accidental background rate was determined by sep- 
arate estimates of the singles rate for both e + -like and 
n-like signals. To minimize possible biases due to the trig- 
ger, both estimates were performed by looking at "iso- 
lated old" events, i.e. at the hits recorded in the NNADC 
event buffer (storing up to 9 events) for which the pre- 
vious and following events were more than 1 ms distant 
(isolated) and which occurred at least 2 ms before the L2 
trigger (old). The event selection was operated by apply- 
ing the same cuts (energy and distance from edge) used 
for the candidate selection; in the 1997 run we found 
R e + = (64.8 ± 0.1) s -1 for the e + -like event rate and 
Rn = (45 ± 2) h" 1 for the n-like one. The rate of acci- 
dental coincidences e + — n in the 2 -f- 100 /is time window 
therefore turns out to be (3.4 ± 0.15) d _1 . An additional 
reduction factor fd = (0.12±0.01) was applied to account 
for the selection operated by the e + — n distance cut. The 
resulting background rate is (0.42 ± 0.05) d _1 , which fits 
with the previous determinations of the total and the cor- 
related background for the 1997 run. 

We checked that the measured single event rates R e + 
and R n are in agreement with estimates based on the ra- 
dioactive contaminants present in the detector and its sur- 
roundings [SUED]. 



9.6 Neutrino yield versus power 

The CHOOZ experiment gave us the unique opportunity 
of measuring the neutrino flux before either reactor started 
working and also after the reactors were shut down. We 
were therefore able to collect enough reactor-off data to 
determine the background rate and could further measure 
the neutrino flux while the reactors were ramping up to 
full power. By fitting the slope of the measured rate ver- 
sus reactor power it is possible to determine the neutrino 
rate at full power, which can then be compared with the 
predicted neutrino production rate at full power. This pro- 
cedure provides not only a very powerful tool to determine 
the neutrino deficit and test the oscillation hypothesis, but 
also gives an independent estimate of the background rate 
to be compared with the values quoted above. 

The fitting procedure is as follows. For each run i the 
predicted number of neutrino candidates Ni is derived 
from the sum of a signal term, linearly dependent on reac- 
tor power Wu and W 2 i, with the background B, assumed 
to be constant and independent of power: 

Ni = {B + W U Y U + W 2i Y 2i )At i , (40) 

where Ati is live time of run i and (Yu, Y 2 i) are the positron 
yields induced by each reactor. These yields still depend 
on the reactor index (even in the absence of neutrino os- 
cillations), because of the different distances, and on run 
number, as a consequence of their different fissile isotope 
composition. It is thus convenient to factor Yki into a func- 
tion Xk (common to both reactors in the no-oscillation 



case) and distance dependent terms, as follows: 

1? 

Y ki = {l+r) ki )-±X k , (41) 
L k 

where fc = 1,2 labels the reactors and the rjki corrections 
contain the dependence of the neutrino counting rate on 
the fissile isotope composition of the reactor core. We are 
thus led to define a cumulative "effective" power according 
to the expression 3 

2 L 2 
W*=J2Wki(l + Vki)j?; (42) 
k=i k 

Ea. (|4l7|) then becomes 

N i = (B + W*X)At i , (43) 

X being the positron yield per unit power averaged over 
the two reactors. 

The burn-up corrections rj ki must be evaluated on a 
run by run basis. Since the GEANT routines would have 
been very time-consuming for such a task, we preferred to 
follow a simpler approach than the procedure explained 
at H9.ll 10000 F e 's per run were generated at each reactor 
according to the respective fuel composition; the kinetic 
positron energy was obtained by the simple relation E u = 
T c + + 1.804 MeV (which comes from © in the limit of 
infinite nucleon mass). The cross section ()10f) was thus 
multiplied by a function S(E U ) to correct for the shift in 
the positron energy scale due to the finite neutron recoil 
effect [10]; this correction can be parametrized as 

8{E V ) = 1 - 0.155 x exp ^ 

We then applied the detector response function (evalu- 
ated by Monte Carlo simulations at several positron ener- 
gies ranging from 0.5 MeV up to 10 MeV) to the positron 
kinetic energy and weighted each event according to the 
positron threshold efficiency. 

We built the likelihood function £ for a set of n runs as 
the joint Poissonian probability of detecting Ni neutrino 
candidates in the i-th run when iVj are expected: 

n 

-ln£= -^\n.P(Ni;Ni) (45) 

i=l 

Searching for the maximum likelihood to determine the 
parameters X and B is then equivalent to minimizing l|45(l . 
The minimization procedure is similar to that used in the 
event reconstruction (see ©. 

Both the average positron yield X and the background 
rate B are assumed to be constant. This is true, by defini- 
tion, for the positron yield, since the effect of the threshold 

3 The "effective" power is the thermal power necessary to 
make the same number of events with a single reactor located 
at the reactor 1 site. It must provide 9.55 GW at full operating 
conditions at the start of reactor operation. 
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resetting on the positron efficiency is already included in 
the rjki correction factors. On the other hand, we find that 
the background rate changes significantly when the trig- 
ger threshold are reset. The complete run sample therefore 
needs to be divided into three periods, according to the 
dates of the threshold resetting, and the fit parameters 
for each period drawn separately. The results are listed in 
Tab.|7| 

Table 7. Summary of the likelihood fit parameters for the 
three data taking periods. 



all data 



period 


f 


2 


3 


starting date 


97/4/7 


97/7/30 


98/1/12 


runs 


579->1074 


1082^1775 


1778->2567 


live time (h) 


1831.3 


2938.8 


3268.4 


reactor-off 


38.9 


539.5 


2737.2 


time (h) 








jWdt (GWh) 


7798 


10636 


2838 


B (O 


1.25 ±0.6 


1.22 ±0.21 


2.2 ±0.14 


X (d^GW- 1 ) 


2.60 ±0.17 


2.60 ±0.09 


2.51 ±0.17 


X'fdof 


136/117 


135/154 


168/184 


n v (d" 1 ) 


24.8 ± 1.6 


24.8 ±0.9 


24.0 ± 1.6 


(at full power) 









By averaging the signal X over the three periods, we 
obtain 



(X) = (2.58 ± 0.07) counts d _1 GW" 



(46) 



corresponding to (24.7 ± 0.7) daily neutrino interactions 
at full power. 

The background rate was stable throughout the 1997 
run, then it suddenly increased by close to a factor of 2 in 
coincidence with the last threshold resetting (see period 3 
column) . This can be explained in terms of a larger acci- 
dental background rate, as a result of threshold lowering. 

Fig. 1491 shows the daily number of neutrino candidates 
as a function of reactor power for all the data sets together. 
The superimposed line corresponds to the average fitted 
signal and background values. 



9.7 Neutrino yield for individual reactors 

The same fitting procedure can be used to determine the 
contribution to the neutrino yield from each reactor and 
for each energy bin of the positron spectra. After splitting 
the signal term into separate yields and introducing a de- 
pendence on the positron energy, Eq. (|43|l can be rewritten 
as 

NiiEj) = (BiE^ + W^X^ + WUE^X^E^AU 

(47) 

The spectrum shape is expected to vary, due to fuel aging, 
throughout the reactor cycle. Burn-up correction factors 
rjki then need to be calculated for each bin of the positron 
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Fig. 49. Daily number of u e candidates, as a function of the 
reactor power. 



spectrum. The fitted yields, averaged over the three pe- 
riods, are listed in Tab. |H] and plotted in Fig. [SD] against 
the expected yield in absence of neutrino oscillations. The 



Table 8. Experimental positron yields for both reactors (Xl 
and X2)and expected spectrum (X) for no oscillation. The 
errors (68% C.L.) and the covariance matrix off-diagonal ele- 
ments are also listed. 



E a + 


Xl±<Tl 


X2 ± a 2 


X 


C12 




MeV 


(counts d^GW" 1 ) 








1.2 


0.151 ±0.031 


0.176 ±0.035 


0.172 


-2.2 10" 


-4 


2.0 


0.490 ± 0.039 


0.510 ±0.047 


0.532 


-1.5 10" 


-4 


2.8 


0.656 ±0.041 


0.610 ±0.049 


0.632 


-3.5 10" 


-4 


3.6 


0.515 ±0.036 


0.528 ± 0.044 


0.530 


-3.3 10" 


-4 


4.4 


0.412 ±0.033 


0.408 ± 0.040 


0.379 


-2.0 10" 


-4 


5.2 


0.248 ±0.030 


0.231 ±0.034 


0.208 


-0.7 10" 


-4 


6.0 


0.102 ±0.023 


0.085 ± 0.026 


0.101 


-1.3 10" 


-4 



yield parameters X\ , Xi are slightly correlated, as shown 
in Tab. [SI this correlation (which does not exceed 20%) is 
always negative, since, at given candidate and background 
rates, an increase in reactor 1 yield corresponds to a de- 
crease in reactor 2 yield (and vice versa). X\ and X-i are 
expected to be the same in the case of no-oscillation. 



9.8 Neutrino yield versus fuel burn-up 

As previously remarked, the contributions of the main fis- 
sile isotopes to the thermal power change in the course of 
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Fig. 50. Positron yields for the two reactors, as compared with 
expected yield for no oscillations. 



an operating period; this should produce a corresponding 
decrease in the total neutrino counting rate as well as a 
modification in the spectral shape. The magnitude of this 
variation amounts to « 10% throughout the first cycle 
of the CHOOZ reactors, as we have already shown (see 
t !2.6(l . which by far exceeds the statistical and systematic 
accuracy of the neutrino flux. We were then forced to fol- 
low the dynamics of the fuel burn-up and to apply daily 
corrections to the thermal power in order to restore the 
linearity with the positron yield (see eas. l|4l"|) . 1(121) ). 

We set out to establish whether the neutrino counting 
rate varies with the reactor burn-up according to predic- 
tions, with the purpose of improving the reliability and 
internal consistency of our results. The number N v of neu- 
trino events, recorded in a time interval At, by a detector 
at a distance L from the core of a reactor working at a 
power W can be derived from Eq. (JJJ ; by inverting we find 



Ef 



WAt N p e 



(48) 



(Tf and Ef being respectively the reaction cross section 
HI Up and the average energy absorbed in the core in a 
single fission. The ratio <jf/Ef, which contains the de- 
pendence on the fuel composition, can be evaluated for 
either reactor by considering the runs where reactors were 
alternatively on; runs are selected if the thermal energy 
is J W dt > 850 GW d (which is the energy released in 
one day by a reactor at 20% of full power). The number 
of neutrino events is obtained by subtracting the num- 
ber of background events (as determined in the previous 
Section) from the number of candidate events collected 
in a run. The resulting crt/Ef values are finally grouped 
in 1000 MWd burn-up intervals and plotted vs. reactor 
burn-up in Fig. I5ir ieft1. 
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Fig. 51. Variation of the measured neutrino counting rate, as 
a function of the fuel burn-up for separate (left) and combined 
(right) reactor data and comparison with predictions. Error 
bars include only statistical uncertainties. 



A more significant test is performed by using combined 
reactor information. Each run is assigned an average burn- 
up defined as 



- s _ PiW l lL\ + p 2 W 2 /Ll 
Wi/Ll + W 2 /L 2 2 



(49) 



where (3\ and [3 2 are the burn-up values for the two reac- 
tors. 

Similarly, aj/Ef is obtained by 



Ef 



= 4vr( 



W x At 



Ll V Ny 

W 2 At> N p e 



(50) 



The resulting values are plotted in Fig. I51f right - ) and com- 
pared with expectations; the compatibility is excellent, 
since \ 2 — 4-69 with 7 d.o.f., corresponding to 69.8% \ 2 ~ 
probability. The data compatibility with a flat distribution 
is also good (x 2 = 5.75 with 6 d.o.f., P x 2 = 45.2%), but 
lower than in the previous comparison. 



9.9 Neutrino direction 

The use of reaction J5J) to detect low energy antineutrinos 
in large volume scintillator detectors is an ideal tool for 
measuring the antincutrino energy spectrum. Moreover, 
it provides a good determination of the antineutrino in- 
coming direction |61ll62j . This determination is based on 
the neutron boost in the forward direction, as a result of 
the kinematics of the reaction; the neutron then retains a 
memory of the source direction, which survives even after 
collisions with the protons in the moderating medium. 

In CHOOZ we exploited this neutron recoil technique 
to locate the reactor direction with the twofold ob- 
jective of testing our event reconstruction method and 
tuning our Monte Carlo simulations of the slowing down 
and capture of neutrons in the scintillator. We then stud- 
ied a possible extension of this technique to much larger 
scintillation detectors, such as Kamland, and their capa- 
bilities in locating astrophysical neutrino sources, such as 
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Supernovae. The neutron emission angle with respect to 
the incident neutrino direction is limited to values below 
~ 55°; this results from Fig. [521 where the neutron an- 
gle 9 n is plotted vs. the neutron kinetic energy T n (which 
extends up to ~ 100 keV). Moreover, the neutron mod- 
eration maintains some memory of the initial neutron di- 
rection [HI]: in each elastic scattering the average cosine 
of the outgoing neutron is (cos# n ) = 2/3A, A being the 
mass number of the scattering nucleus. The direction is 
thus best preserved by collisions on protons, which is also 
the most effective target nucleus at energies below 1 MeV. 
The neutron mean free path rapidly reduces during mod- 
eration, since the scattering cross section rapidly increases 
at increasingly lower neutron energies; so the bulk of the 
neutron displacement is due to the first two or three colli- 
sions. The isotropic diffusion of thermalized neutrons does 
not affect the average neutron displacement along the neu- 
trino direction. 



Table 9. Measurement of neutrino direction: data and Monte 
Carlo. The angle 8 is the la— uncertainty on 77 c direction. 




NeutroiLKinecc Energy- (MeV) 



Fig. 52. Neutron emission angle (with respect to the incident 
V e direction) vs. its kinetic energy; the discrete structure of 
lower-left part of the picture is an effect of the logarithmic 
scale for the abscissa combined with the V a energy binning. 



The average neutron displacement in CHOOZ has been 
calculated to be 1.7 cm. Since the experimental position 
resolution is a x w 19 cm for the neutron and the collected 
neutrino statistic is « 2500, the precision of the method is 
w 0.4 cm; so neutron displacement can be observed at ~ 
4cr level. The average direction (in spherical coordinates) 
of the two reactors in the CHOOZ detector frame was 
measured by standard surveying techniques to be cf> = 
(-50.3 ± 0.5)° and 9 = (91.5 ± 0.5)°. 

The neutrino direction resulting from the data and the 
associated uncertainty, is presented in Tab. |UJ the Monte 
Carlo predictions, for a sample with the same statistics, 
is also listed for comparison. The measured direction has 
a 16% probability of being compatible with the expected 
value, while the probability of a fluctuation of an isotropic 
distribution is negligible. The average neutron displace- 
ment was found to be (1.9 ± 0.4) cm, in agreement with 
what expected. 
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10 Neutrino oscillation tests 

Three different approaches were adopted to analyse the 
results of the experiment. The first one ("analysis A" in 
what follows) uses the predicted positron spectrum in ad- 
dition to the measured spectra for each reactor. This ap- 
proach uses the absolute normalization and is therefore 
sensitive to the mixing angle even for large Sm 2 values, 
where the oscillation structure could no longer be resolved 
in the energy spectrum and the oscillation limits exclu- 
sively depend on the uncertainty in the absolute normal- 
ization. 

The second approach ("analysis B") uniquely relies 
on a comparison of measurements taken at different re- 
actor core-detector distances. We are led to limits on the 
oscillation parameters which are practically independent 
of the uncertainties in the reactor antineutrino flux and 
spectrum; other major sources of systematic uncertain- 
ties, such as detection efficiencies and the reaction cross 
section, also cancel out. The result of this analysis can 
thus be regarded as free from all these systematic uncer- 
tainties. 

The third approach ( "analysis C" ) is somewhat inter- 
mediate between the first two analyses. It uses the shape 
of the predicted positron spectrum, while leaving the ab- 
solute normalization free. The only contribution to the 
systematic shape uncertainty comes from the precision of 
the neutrino spectrum extraction method 26 . 



10.1 Global test (Analysis A) 

In a simple two-neutrino oscillation model, the expected 
positron yield for the k-th reactor and the j-th energy 
spectrum bin, can be parametrized as follows: 



X(E. 3 ,L k ,9,bm 2 ) = X(E 3 )P(E 3 ,i, k 



1, 



L k ,9,Sm 2 ), 
,7 k = l,2) 



(51) 



where X(Ej) is the distance- independent positron yield in 
the absence of neutrino oscillations defined in the previ- 
ous section, L k is the reactor-detector distance and the 
last factor represents the survival probability averaged 
over the energy bin and the finite detector and reactor 
core sizes. The procedure to compute such a probability 
at varying oscillation parameters is similar to the method 
used to calculate the burn-up corrections to the positron 
yields (see H9.6f) : the positron spectrum is obtained by 
Ea.lpT7)l after adding the detector response function and 
the size function; the same procedure was applied to ob- 
tain the spectrum for no oscillations; the probability P(Ej , L k 



, Sm 
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results then from the ratio of the j-th bin contents of the 
two spectra. 

In order to test the compatibility of a certain oscilla- 
tion hypothesis (9,8m 2 ) with the measurements, we built 
a x 2 statistic containing the experimental yields for the 7 
energy bins at the two positions Lk are listed in Tab.|Hl We 
grouped these values into a 14-element array X arranged 
as follows: 



X=(X 1 (E 1 ),...,X 1 (E r ),X 2 (E 1 ) 



,X 2 (E 7 )), (52) 



and similarly for the associated variances. These compo- 
nents are not independent, as already noted in H9.7I By 
combining the statistical variances with the systematic un- 
certainties related to the neutrino spectrum, the 14 x 14 
covariance matrix can be written in a compact form as 
follows: 



V t3 = 5 i . J (a l + a i ) + (Sij-r + 8 i>j+r )a\ 



(i,j = l,...,14), (53) 

where cr, are the statistical errors associated with the yield 
array (|52l) , <jj are the systematic uncertainties and are 
the covariance of reactor 1 and 2 yield contributions to 
the i-th energy bin (see Tab. [HJ- These systematic errors, 
including the statistical error on the measured /3-spectra 
measured at ILL 26 as well as the bin-to-bin system- 
atic error inherent in the conversion procedure, range from 
1.4% at 2MeV (positron energy) to 7.3% at 6MeV and 
are assumed to be uncorrelated 4 . 

We still have to take into account the systematic error 
related to the absolute normalization; combining all the 
contributions listed in Tab. 1101 we obtain an overall nor- 
malization uncertainty of a a — 2.7%. We may define the 
following x 2 function 

X 2 (9,8m 2 ,a,g) = 



i=l j=l 

Vr j 1 (x 3 -aX{gE ,L ] ,e,6m 

2 



a — 1 



.9-1 



(54) 



where a is the absolute normalization constant, g is the 
energy-scale calibration factor, Lij = L\ for i,j<7 and 
Lij = L 2 for i,j > 7. The uncertainty on g is a g = 1.1%, 
resulting from the accuracy on the energy scale calibra- 
tion (16keV at the 2.11 MeV visible energy line associ- 
ated with the n-capture on Hydrogen) and the 0.8% drift 

4 The extraction of the neutrino spectra from the (5 mea- 
surement at ILL should introduce a slight correlation of the 
bin-to-bin error systematic error. Nevertheless, the overall un- 
certainties on the positron yields are dominated by statisti- 
cal errors, so that neglecting the off-diagonal systematic error 
matrix does not affect the oscillation test significantly. Also 
previous reactor experiments, even with much lower statistical 
errors, did not take this correlation into account. 



in the Gd-capture line, as measured throughout the ac- 
quisition period with high-energy spallation neutrons (see 
Fig. I38fl . The x 2 m jUt thus contains 14 experimental 



Table 10. Contributions to the overall systematic uncertainty 
on the absolute normalization factor. 



parameter 



relative error (%) 



reaction cross section 


1.9% 


number of protons 


0.8% 


detection efficiency 


1.5% 


reactor power 


0.7% 


energy released per fission 


0.6% 


combined 


2.7% 



errors with 2 additional parameters, yielding altogether 
16 variances. The x 2 value for a certain parameter set 
(9,5m 2 ) is determined by minimizing Ij54(l with respect to 
the gain factor g and to the normalization a; the mini- 
mization then leads to 12 degrees of freedom. The mini- 
mum value Xmin = 5.0, corresponding to a % 2 probability 
P x = 96%, is found for the parameters sin 2 (29) = 0.23, 
5m 2 = 8.1 ■ 10" 4 eV 2 , a = 1.012, g = 1.006; the result- 
ing positron yields are shown by solid lines in Fig. 1531 and 
superimposed on the data. Also the no-oscillation hypoth- 
esis, with x 2 (0, 0) = 5.5, a — 1.008 and g — 1.011, is found 
to perfectly fit with the data (P x = 93%). 
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Fig. 53. Positron yields for reactor 1 and 2; the solid curves 
represent the predicted positron yields corresponding to the 
best-fit parameters, the dashed one to the predicted yield for 
the case of no oscillations. 



To test a particular oscillation hypothesis (9, Sm 2 ) 
against the parameters of the best fit and to determine the 



:S(> 



90% confidence belt, we adopted the Feldman & Cousins 
prescription |65| . The "ordering" principle is based on the 
logarithm of the ratio of the likelihood functions for the 
two cases: 



by noting that 



X(9,Sm 2 ) = X 2 (9,Sm 2 ) 



v 2 . 

Amm 



(55) 



where the mimimum \ 2 value must be searched for within 
the physical domain (0 < sin 2 (29) < 1, Sm 2 > 0). Smaller 
A values imply better compatibility of the hypothesis with 
the data. The A distribution for the given parameter set 
was evaluated by performing a Monte Carlo simulation of 
a large number (5000) of experimental positron spectra 
whose values are scattered around the predicted positron 
yields X(Ei, Li, 9, Sm 2 ) with Gaussian- assumed variances 
<7j and correlation coefficients given by (|53|l . For each set 
we extracted the quantity X c (9, Sm 2 ) such that 90% of the 
simulated experiments have A < A c . The 90% confidence 
domain then includes all points in the (sin 2 (29), Sm 2 ) plane 
such that 



X exp (9,Sm 2 ) < X c (9,Sm 2 ), 



(56) 



R ps [l 



1.27<5m 2 ii, 

■ El ' 



.1.275m 2 L 2 , 2l 



l-2( 



E v 
1.27 'Sm 2 , 



"LSL, (57) 



where L is the average reactor-detector distance and SL 
is the difference of the two distances. Therefore an exper- 
iment which measures this ratio with an uncertainty a is 
sensitive to oscillations down to mass-difference values as 
low as 

Ev " (58) 



Sm 



1.27 V 2LSV 



k being the number of standard deviations corresponding 
to the chosen confidence level. This value can be com- 
pared to the sensitivity limit Sm^ inherent in analysis A 



Sm 2 



2lL 6m ° 



2Sm% 



1.5- lO^eY^ 



(59) 



Although twice as large, this limit is lower than the lowest 
Sm 2 value allowed by Kamiokande (see Fig. l55|) . The ratio 
R(Ei) = Xi(Ei) I 'X 2 (Ei) of the measured positron yields 
must be compared with the expected values; since the ex- 
pected yields are the same for both reactors in the case of 
no-oscillations, the expected ratio for the i-th energy bin 
reduces to the ratio of the average survical probability in 
that bin. We can then build the following \ 2 statistic: 



/ R(Ei) — R(E t ,9, Sm 2 ] 



SR(Ei) 



(60) 



where X exp is evaluated for the experimental data for each 
point in the physical domain. 

The acceptance domain at the 90% C.L. (solid line) 
and 95% C.L. are shown in Fig. ESI all the parameters ly- 
ing to the right of the curves are excluded by CHOOZ with 
the indicated confidence level, while the parameter regions 
on the left are still compatible with our data. The region 
allowed by Kamiokande for the — > v c oscillations is also 
shown for comparison; this hypothesis, a possible expla- 
nation for the deficit in the atmospheric neutrino flux, 
is excluded. The Sm 2 limit at full mixing is 7 • 10 _4 eV 2 ; 
the limit for the mixing angle in the asymptotic range of 
large mass differences is sin 2 (29) = 0.10. 



10.2 Two-distance test (analysis B) 

The predicted ratio of the two-reactor positron yields equals 
the ratio of the corresponding survival probabilities. At 
full mixing (sin 2 (20) = 1) and at low mass differences 
(Sm 2 w 10~ 3 eV 2 ), this ratio can be approximated by 



where SR(Ei) is the statistical uncertainty on the mea- 
sured ratio. The minimum x 2 value (Xmin — 0-78 over 5 
d.o.f.) occurs at sin 2 (26 | ) = 1 and Sm 2 — 0.6 eV 2 ; the com- 
patibility of the no-oscillation hypothesis is still excellent 
(see Fig. EU, as x 2 (0, 0) = 1.29. 
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Fig. 54. Measured ratio of experimental positron yield, com- 
pared with the predicted ratio in the best oscillation hypothesis 
(solid line) and in the case of no oscillations (dashed line). 



We adopted the same procedure described in the pre- 
vious section to determine the confidence domain in the 
(sin 2 (29), Sm 2 ) plane and for each point in this plane we 
simulated the results of 5000 experiments. If the positron 
yields of both reactors are Gaussian distributed around 
the predicted values X i:2 , and if X 2 /<r 2 is sufficiently large 
so that the X 2 generated have only positive values, then 
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the variable 



Z 



Xi 

x 1 



X2 
x 2 



X\ — X2R 




alR 2 - 2a 12 R 



(61) 



is normally distributed with zero mean and unit vari- 
ance |56j . We could then use a normal random generator 
to extract Z and invert (|61|) to determine the ratio R(Ei). 

The exclusion plot obtained at 90% C.L. is also shown 
in Fig. 1551 Although less powerful than the previous anal- 
ysis, the region excluded by this oscillation test almost 
completely covers the one allowed by Kamiokande. 



10.3 Shape test (Analysis C) 

This test is similar to analysis A, the only difference be- 
ing related to the hypothesis on the absolute normaliza- 
tion. In analysis A we fixed the integral counting rate to 
be distributed around the predicted value (a = 1), with 
<j a — 2.7% systematic uncertainty; in the shape test, on 
the contrary, we gave up any constraint on the normaliza- 
tion parameter (which is equivalent to having a a = 00). 
The x 2 statistic for this test is the same expression as (|54[1 
without the term depending on the normalization, so we 
can write 

X 2 {0,6m 2 ,a,g) = 

14 14 



^2^(Xi- aX(gE h Li,0,6r 
i=i j=i 

V£ 1 (Xj - aX(gE j ,^,6, 5m 2 



9-1 



(62) 



This x 2 has a minimum value Xmin = 2-64 (over 11 de- 
grees of freedom) at sin 2 (26>) = 0.23, 6m 2 = 2.4 • 10~ 2 eV 2 
and g = 1.008; the null hypothesis gives instead x 2 (0, 0) = 
5.5 with g — 1.006. The exclusion plot, obtained accord- 
ing to the Feldman-Cousins prescriptions, is shown in 
Fig. [SB] and compared to the results of the other tests. 
As in the case of the two-distance test it is not sensi- 
tive to oscillations at large squared mass-difference values 
(5m 2 ^O.leV 2 ), where the oscillation length 



L osc (m) 



2.48E^(MeV) 
8m 2 (eV 2 ) 



100 m 



(63) 



becomes much lower than the average reactor-detector 
distance. The sin 2 (29) limit at the maximum oscillation 
probability is similar to that obtained with analysis A. 



10.4 Implications of the CHOOZ results 

The importance of the CHOOZ results on neutrino oscil- 
lations has been pointed out by many authors 0©,66 




analysis A 
analysis B 
analysis C 
1% CL Kamiokande (multi-GeV) 
>% CL Kamiokande (sub+multi-GeV) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin 2 (28) 



Fig. 55. Exclusion plot contours at 90% C.L. obtained by the 
three analyses presented. Analysis A uses both shape and nor- 
malization of the background-subtracted positron spectrum; 
analysis B uses the baseline difference between the two reac- 
tors; analysis C uses the spectrum shape. The Kamiokande 
t>n «-» v e allowed region to atmospheric neutrino oscillations is 
also shown for comparison. 



]. In a three- flavour neutrino mixing frame, 
the see-saw mechanism indicates a mass hierarchy (m\ 
777-2 <C 7773) from which 



(S777 2 2 <C Sm 2 3 . 



(64) 



The only available possibile explanation for solar and 
atmospheric neutrino anomalies through neutrino oscil- 
lations is that Smi2 is relevant for the transitions of solar 
neutrinos and 8m\z is the one probed by atmospheric and 
L-B accelerator neutrino experiments. By approximation 
(|64|l it is possible to show that the CP-violating phase in 
the 3x3 mixing matrix does not give rise to observable 
effects, and that the mixing angle #12 (associated with the 
lower mass states) can be rotated away in the atmospheric 
neutrino analysis. Therefore, condition (|64|l implies that 
the transition probability in atmospheric and L-B experi- 
ments depends only on the largest mass squared difference 
<5tt7,23 an d the elements U a z connecting flavour neutrinos 
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Fig. 56. Results of a three flavour mixing analysis of separate and combined Super-Kamiokande and CHOOZ data, for five 
representative values of m 2 = 8m 2 23- The analysis concerns the 79.5 kTy data sample for Super-Kamiokande and the first 
CHOOZ result. 



with v 3 \ it assumes then the simple form 
,(1 



P — < 



1 - 
if a 



4E& v 



U 2 a3 ) 



l.27Sm 2 {eV 2 )L(m) 
E(MeV) 



[if a ^ 0. 



1.27(5m 2 (eV 2 )L(m) 
E{ MeV) 



(65) 



In particular, the survival probability reduces to the usual 
two-flavour formula with sin 2 26> = 4C/ 2 3 (1 — U^ 3 ). 
Therefore, information on the parameter U e3 can be ob- 
tained from the CHOOZ exclusion plot. The upper limit 
for sin 2 29 implied by the exclusion plot in Fig.[!)3]is sin 2 29 
< 0.1 for 5m 13 > 2 • 1CT 3 ; it follows that 



C/ 2 3 < 0.03 



U 2 3 > 0.97 



(66) 
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Large values of values of U e 3 , those allowed by the second 
inequality in Ea. (|66|) . are excluded by the solar neutrino 
data pP; the solar v c survival probability would be larger 
than 0.95, which is incompatible with the deficit observed 
in all solar neutrino experiments. The CHOOZ result then 
constrains the mixing of electron neutrinos with all other 
flavours, in the atmospheric and L-B range, to small val- 
ues. 

It has been remarked that even the most recent 
Super-Kamiokande data alone do not completely exclude 
sizable v c mixing. This can be seen in Fig. 1561 where 
the confidence regions for separate and combined Super- 
Kamiokande and CHOOZ data are shown for different 
values of m 2 = <5m 2 23 (the mass difference relevant for 
atmospheric neutrinos). Each set of mixing parameters 
(U e 3,Ufj,3,U T 3) is associated with a point embedded in 
a triangle graph (equilateral due to the unitarity condi- 
tion), whose corners represent the flavour eigenstates; by 
convention v is assigned the upper corner. The distance 
from the three sides are equal to U 2 3 , U 2 3 , U 2 3 ; in a two- 
flavour scheme (£/ a3 = 0), the point is bound on the side 
connecting the two mixed flavour states (and the mean 
point of that side is associated with the maximum mixing 
hypothesis) . 

The CHOOZ result excludes large horizontal stripes 
in the triangle plot, according to Ea. lffIB*)) : the stripe be- 
comes increasingly narrower at lower Sm 2 23 values, as 
a consequence of the reduced sensitivity to sin 2 29. The 
Super-Kamiokande allowed region lies on the triangle base 
and protrudes towards the centre, which implies a non- 
negligible — > v c oscillation probablility. Yet in the com- 
bined analysis graph the allowed region is significantly 
flattened on the base, thus indicating a dominance of the 
v ^ — > v T maximum mixing hypothesis. 

The impact of CHOOZ in constraining the v c mixing is 
more evident in the bilogarthmic plot of Fie . 1571 where the 
allowed regions (90 and 99 % C.L.) in the m 2 (= <5ra 2 2 3), 
tan 2 (j) (<p = 613) parameter space are shown for the Super- 
Kamiokande data only (less restrictive limits) and for the 
combined Super-Kamiokande and CHOOZ data. The im- 
provement obtained by using the CHOOZ data is of about 
one order of magnitude. 



11 A novel method for the derivation of 
confidence regions 

The derivation of confidence regions for parameters, in 
neutrino oscillation experiments, has been the focus of 
attention and debate in recent years. There has been a 
general consensus in applying the Feldman and Cousins 
method (FC), which, apart from its intrinsic merit, allows 
a simpler comparison between experiments. Discussions 
about the statistical analysis are far from conclusive. We 
therefore tested an interesting new scheme which may rep- 
resent a further step towards an objection-free procedure. 
It tends to produce confidence regions which are some- 
what larger in size and, similarly, its application would 
more or less affect the results of all experiments. The new 
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Fig. 57. 90% (solid line) and 99% C.L. (dotted line) limits in 
the m 2 (= <5m 2 23), tan 2 (f> (</> = #13) parameter space obtained 
with the use of Super-Kamiokande alone (less restrictive lim- 
its) and CHOOZ +Super-Kamiokande. 



prescription |71j in confidence interval estimate is based 
on the concept of strong Confidence Level (sCL) , to make 
an inference on neutrino oscillation parameters. This new 
method is purely frequentist, just like the FC "unified 
approach" , yet it seems to be free from any pathologies 
inherent in every frequentist procedure proposed to date. 
Moreover, the definition of strong limits can be generalized 
in a natural way to include exact frequentist treatment of 
systematics [72]. Previous limits (Fig. I55f) were computed 
by replacing the likelihood function with the one obtained 
by maximizing with respect to the systematic parame- 
ters ("profile" likelihood). This is an approximation |75] 
often used to overcome computational difficulties. It has 
been shown [72] that this approximation no longer holds 
in cases where the systematic error is of the same order of 
magnitude of the statistical one. 

11.1 Application of the method 

In this method we need to construct the probability den- 
sity function (pdf) of our experiment in terms of the data 
(energy spectra) and the parameters (oscillation + nor- 
malization). For practical reasons, we decided to condense 
all the information (shape + normalization) contained in 

a spectrum X into the three parameters (sin 2 (2#), 8m 2 ), 



40 



a) for which the likelihood is maximum; any parameter 
set uniquely identifies a particular energy spectrum. The 
set of maximum likelihood estimators (MLE) is our ob- 
servable. 

In order to obtain strong confidence limits, we com- 
puted the probability density function numerically. The 
domain of oscillation parameters (0 < sin 2 (26*) < 1, 10 -4 
< Sm 2 < leV 2 ) was sampled by a 100 x 100 grid (with 
a constant binning in log 8m 2 ); the range 0.95 < a < 
1.05 was divided into 10 cells 5 . For each set of (sin 2 (26*), 
5m 2 , a) we generated 5 • 10 4 spectra scattered around 
aX (sin 2 (28) , 5m 2 ) with uncertainties given by the covari- 
ance matrix. MLE for each experiment are searched for 
in the same space spanned by the parameters but with a 
coarser grid (20 x 20 x 10). 
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Fig. 58. Exclusion plot at 90% sCL for the oscillation parame- 
ters based on the differential energy spectrum; the FC contour, 
obtained with "correct systematics" treatment, is also shown. 



12 Conclusions 

The CHOOZ experiment stopped taking data in July 1998, 
about 5 years after the submission of the proposal for ap- 
proval. With more than 1-year data taking, the statistical 
error (2.8%) on the neutrino flux matched the goal (3%) 
of the proposal. Accurate estimatea of the detection effi- 
ciencies as well as precise measurements of the detector 
parameters also allowed us to keep the systematic uncer- 
tainty (2.7%) below expectations (3.2%). 

We found (at 90% confidence level) no evidence for 
neutrino oscillations in the V e disappearance mode, for 
the parameter region given by approximately 5m 2 > 7 ■ 
10 _4 eV 2 for maximum mixing, and sin 2 29 = 0.10 for 
large 8m 2 . Less sensitive results, based only on the com- 
parison of the positron spectra from the two different- 
distance nuclear reactors (and therefore independent of 
the absolute normalization of the V e flux, the number of 
protons and the detector efficiencies) were also presented. 

Our result does not allow the atmospheric neutrino a- 
nomaly to be explained in terms of — > v c oscillations, 
thus leaving, in a three-flavour mixing scheme, the — > 
v T possibility. 

Many cross-checks were performed on the data to test 
the internal consistence and improve the reliability of our 
results. As a by-product, we have shown that the use of 
reaction allowed us to locate the V c source within a 
cone of half-aperture ~ 18° at 68% confidence level. 
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scintillator Hydrogen content. 



The results of our computation are shown in Fig. 1581 
The confidence bounds obtained are significantly higher 
than those obtained by the procedure explained in the 
previous section. In fact, oscillations v c — > 7/ x are ex- 
cluded for 5m 2 > 8 • 10" 4 eV 2 at maximum mixing and 
sin 2 (20) > 0.17 at large 5m 2 values. It should be noted 
that the limits quoted are only slightly looser than those 
obtained by using the FC prescription with the correct 
inclusion of systematics, as shown in Fig. 1581 



5 We verified that values of a outside the considered range 
give no further contribution to the projection of the band on 
(sin 2 (2#), <5m 2 )-space. 
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